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PREFACE 


This  report  was  prepared  by  M.  K.  Wahi,  S.  M.  Warren,  and  H.  H.  Straub  of  the  Boeing 
Commercial  Airplane  Company  under  a USAF  Contract  F33657-74-C-0J  29  (extended  ).  The 
program  was  divided  into  three  tasks.  Task  I involved  the  formulation  of  a tire  correlation 
model,  a test  outline  to  validate  the  model,  establishment  of  a friction  prediction  subsystem 
specification  criteria  and  evaluation  of  existing  ground  vehicles.  The  work  completed  under 
Task  1 of  this  contract  was  performed  from  May  1975  to  December  1975  and  all  aspects  of 
that  work  are  described  in  a separate  report  ASD-TR-77-7.  .. , . 

Task  !i  involved  a sensitivity  analysis  of  airplane  braking  distance  by  using  the  Boeing  Brake 
Control  Simulator  for  the  USAF  B-52,  KC-135  and  F-l  1 1 airplanes  to  validate  the  general 
prediction  model  developed  in  the  previously  contracted  effort.  Task  111  objectives  were  to 
establish  compatibility  between  Task  11  and  Task  1 subsystems  and  to  recommend  a test 
program  to  verify  the  effectiveness  and  reliability  of  the  Total  Braking  Prediction  System 
(TBPS).  The  present  volume  describes  all  essential  aspects  of  the  work  performed  in  com- 
pleting Tasks  11  and  111  of  this  contract.  Volume  11  describes  the  hardware  and  antiskid 
systems  used  on  the  brake  control  simulator  as  well  as  the  test  conditions  and  parameters 
used  in  developing  data  required  for  the  dimensional  analysis.  The  work  described  herein 
was  performed  from  August  1975  to  December  1976. 

The  authors  are  indebted  to  Messers  N.  S.  Attri.  and  A.  J.  P.  Lloyd  for  their  guidance  and 
technical  contribution  as  respective  program  managers  at  various  stages  of  the  contract. 

T.is  authors  are  also  indebted  to  Messers  W.  B.  Tracy  ASD/ENFEM.  D.  B.  Tremblay  ASD./ 
ENFE.  Lt.  Coi.  R.  Kennah.  ASD/AEAA  and  Maj.  R.Cauley,  AFCEC/EMR  of  the  USAF  for 
their  program  support 
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SUMMARY 


An  airplane  braking  sensitivity  study  has  been  extended  to  the  F-l  1 1A,  KC-135  and  B-52H 
airplanes  with  an  analog-hardware  brake  cortrol  simulation.  As  indicated  in  earlier  studies, 
the  peak  available  ground  friction  (p),  drag  device  effectiveness  brake  application 

speed  (V).  air  density  (p).  and  engine  idle  thrust  (Fe)  are  the  resulting  parameters  which 
have  the  most  significance  on  braking  distance.  The  braking  distance  prediction  equations 
involving  these  parameters  have  been  developed  for  each  of  the  above  three  airplanes,  and 
are  compatible  with  those  developed  for  other  aircraft.  Recommendations  for  the  opera- 
tional application  of  the  Total  Braking  Prediction  System  are  described  and  involve  the  use 
of  a suitable  ground  friction  vehicle  (Friction  Prediction  Subsystem),  a proven  tire  correlation 
model,  and  .In  prediction  equations.  The  feasibility  of  the  Total  Braking  Prediction  System 
must  be  verified  with  a recommended  test  program. 


SECTION  I 
INTRODUCTION 


This  report  presents  in  two  volumes  the  results  of  Task  1!  and  III  of  Combat  Traction  II, 
Phase  II,  Extended.  Results  of  Task  1 of  this  program  are  presented  in  a separate  report- 
“Tire  Runway  Interface  Friction  Prediction  System”  ASD-TR-77-7.  The  Task,  summarized 
diagrammatically  in  Figure  1 , describes  a tire  correlation  model  concept,  recommends  a 
test  program  for  validation  of  the  correlation  model,  and  presents  a specification  for  a 
ground  friction  measuring  vehicle. 

Task  II  consists  of  an  airplane  braking  sensitivity  study  similar  to  that  of  References  ! and  2 
as  applied  to  the  F-l  1 1 A,  KC-135  and  B-52H.  Results  of  the  sensitivity  analysis  (Task  JIa, 
Figure  2)  combined  witit  those  of  the  earlier  sensitivity  study  (Reference  1 ) were  used  to 
form  a composite  group  of  significant  factors  affecting  braking  performance..  This  list  of 
factors  was  used  to  form  an  expanded  Braking  Performance  Prediction  Model  (Task  Ub). 
Using  this  expanded  model  individual  prediction  equations  for  each  aircraft  under  consider- 
ation have  been  developed  (Task  lie). 


Results  of  the  tire-runway  friction  study  (Task  I)  and  the  airplane  braking  study  (Task  II) 
were  combined  in  Task  III  to  form  a Total  Airplane  Braking  Prediction  System.  In  addition 
a test  program  is  lecommended  to  verify  the  total  airplane  braking  prediction  system  concept. 


Overall  Program  Plan 
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SECTION  II 

TASK  II  SENSITIVITY  STUDY 


The  objective  of  this  sensitivity  study  was  to  determine  the  parameters  that  influence  the 
broking  distance  of  the  F-l  1 1 A.  KC-135  and  B-52H.  Only  those  parameters  which  vary 
during  the  normal  operation  ot  the  airplane  were  considered.  Parameter  variat.t>!is  involving 
modification  of  the  brake  system  or  dependent  upon  pilot  input  were  not  considered. 

The  data  used  in  the  brake  control  simulator  for  the  three  aircraft  were  supplied  primarily 
by  the  USAF.  A list  of  documents  used  during  data  acquisition  may  be  found  in  the  refer- 
ence list  (References  3 thru  23).  The  data  in  many  cases  were  reduced  to  the  form  required 
for  use  in  the  simulator. 

As  a starting  point,  a baseline  airplane  was  defined  for  each  aircraft.  The  baseline  airplane 
represents  an  aircraft  in  a three-point  taxi  attitude  and  of  typical  (mid-range)  landing  weight, 
approach  speed,  center-of-gravity  location,  landing  flap  setting,  and  engine  thrust.  The  actual 
parameters  required  for  the  airplane  simulator  are  defined  in  Table  1 . Table  2 lists  the  base- 
line values  used  for  each  airplane.  In  some  cases  insufficient  data  was  available  to  define  an 
airplane  parameter,  in  these  cases  a value  was  assumed  based  on  experience  and  engineering 
judgment.  Table  3 lists  the  parameters  (for  each  aircraft)  for  which  an  assumption  was  nec- 
essary. As  pointed  out  in  the  table,  brake  parameters  and  strut  damping  values  were  most 
often  undefined. 

During  the  sensitivity  study  each  parameter  was  changed  and  the  new  value  of  braking  d' stance 
was  evaluated  on  the  simulator.  Some  of  the  variables  were  not  independent,  and  those  groups 
of  interrelated  parameters  were  varied  appropriately  together.  An  example  of  this  is  stall 
speed  and  gross  weight.  During  previous  work  the  range  over  which  a parameter  was  varied 
represented  a normal  operating  range,  however  for  this  study  parameters  were  varied  to  their 
maximum  (minimum)  allowable  operating  values. 
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Table  1 -Definition  of  Simulator  Parameters 


1 PARAMETER  I 

DEFINITION  ”1 

UNIT 

AW 

Effective  wing  area 

CD 

Drag  coefficient 

COD 

Drag  coefficient  without  spoilers 

CL 

Lift  coefficient 

CLD 

Lift  coefficient  without  spoilers 

FEO 

Engine  idle  Thrust  at  zero  velocity 

lbf 

MB 

Height  of  CG  above  ground 

ft 

IYY 

Mass  moment  of  inertia,  pitch 

ft-lb-sec* 

KE 

Change  of  idle  thrust  with  velocity 

lbf-sec/ft 

LA 

Nose  gear  to  CG  distance 

ft 

LB 

Main  gear  to  CG  distance 

ft 

NB 

Number  of  brakes  per  main  strut 

NBA 

Number  of  main  gear  brakes  per  airplane 

NBN 

Number  of  nose  gear  wheels 

NS 

Number  of  main  gear  struts  per  airplane 

RHO 

Air  density 

lbf-sec/ft^ 

VI 

Initial  airplane  velocity 

ft/sec 

VSTOP 

Final  airplane  velocity 

ft/ sec 

UA 

Weight  of  airplane 

lbf 

1 BRAKE  PARAMETERS  1 

KP 

Torque  peaking  gain 

MB 

Mass  of  brake  heat  sink 

lbm 

OMGP 

Wheel  velocity  at  start  of  torque  peaking 

rad/ sec 

PC 

Retractor  spring  pressure 

psi 

TBG 

Torque  gain 

ft-lbf/psi 

THBB 

Temperature  at  initiation  of  fade 

°F 

WN 

Natural  frequency  of  torque  response 

Hz 

ZETA 

Damping  ratio  of  torque  response 

| TIRE  PARAMETERS  ] 

1 0 

Tire  diameter 

in. 

Tire  deflection 

in. 

!M 

Mass  moment  of  inertia  of  tire,  wheel  and  brake 

ft-lbf-sec2 

4 
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Table  1 .—Definition  of  Simulator  Parameters  (Concluded) 


PARAMETER 

DEFINITION 

UNIT 

j TIRE  PARAMETERS  (Continued) 

PI 

Tire  operating  Inflation  pressure 

psi 

PR 

Tire  rated  inflation  pressure 

psi 

RR 

Tire  rolling  radius 

ft 

RT 

Tire  torque  radius 

ft 

| STRUT  PARAMETERS 

CO 

Main  gear  vertical  damping  coefficient 

lbf -sec/ft 

CON 

Nose  gear  vertical  damping  coefficient 

lbf-sec/ft 

CS 

Main  gear  fore-aft  damping  coefficient 

Ibf-sec/ft 

CT 

Torsional  Damping  between  strut  and  brake 

Ibf-ft-sec 

IS 

Mass  moment  of  inertial  of  main  gear  strut 

ft-lbf-sec2 

KO 

Main  gear  vertical  stiffness 

Ibf/ft 

KON 

Nose  Gear  vertical  stiffness 

Ibf/ft 

KS 

Main  gear  fore-aft  stiffness 

lbf/ft 

t 

Effective  strut  length 

ft 

MS 

Effective  mass  of  strut 

Ibf-sec2/ft 

Table  2— Baseline  Values  Used  in  Airplane  Simulation 


Table  3.— Assumed  Airplane  Parameters 


B-52 

KC-I35 

F-111A 

AIRPLANE  PARAMETERS 

a 

CD 

BRAKE  PARAMETERS 

KP 

KP 

KP 

OMGP 

OMGP 

OMGP 

THBB 

THBB 

THBB 

UN 

UN 

UN 

ZETA 

ZETA 

ZETA 

STRUT  PARAMETERS 

CO 

CO 

CO 

CON 

CON 

CON 

CS 

CS 

CS 

CT 

CT 

CT 

SECTION  III 

SENSITIVITY  STUDY  TEST  CONDITIONS 


The  sensitivity  study  involved  changing  a parameter  or  group  of  parameters  to  analyze  the 
resulting  influence  on  braking  distance.  The  test  conditions,  as  listed  in  Table  4,  were  per- 
formed on  the  analog-hardware  simulator  to  determine  braking  distance.  Table  5 lists  the 
parameters  changed  and  the  numerical  value  of  the  parameters  associated  with  each  test 
condition  for  the  three  aircraft. 

Nine  sensitivity  tests  were  formulated  to  analyze  the  sensitivity  of  the  various  aircraft  systems 
to  a parameter  change.  The  tests  were  performed  selectively  at  each  test  condition.  I he  nine 
sensitivity  tests  liave  been  divided  into  three  major  categories.  The  categories  and  the  tests 
associated  with  each  are  as  follows: 

• Stability  studies. 

Test  1 - strut  stability 

• Pei formance  studies: 

Test  2 - touchdown  dynamics 
Test  3 - stabilized  landing 
Test  4 - mu  steps 
Test  5 - wet  runway 

• Hydraulic  system  studies: 

Test  , - frequency  response 
Test  7 - step  response 

Test  8 - antiskid  valve  characteristics 

Test  9 - brake  pressure  - volume  characteristics 

A detailed  description  of  the  test  procedure  and  sequence  can  be  found  in  Volume  11.  The 
three  major  categories  of  sensitivity  tests  are  briefly  described  below  to  point  out  their  general 
significance. 


1.  STABILITY  STUDIES 

System  stability  is  directly  related  to  stopping  performance.  Severe  instability  can  result  in 
the  loss  of  braking  and  can  cause  serious  safety  hazards.  The  study’s  purpose  was  to  evaluate 
the  ability  of  a brake  control  system  to  contribute  to  the  fore-aft  stability  of  the  gear. 
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Table  4.— Test  Conditions 


Baseline  Study 

Nominal  values  of  all  parameters 
Parametric  Studies 
Airplane 

1.  Weight 

a.  Maximum,  with  VI  effect 

b.  High  intermediate  with  VI 

c.  Low  intermediate  with  VI 

d.  Minimum  with  VI 

e.  Maximum  without  VI  effect 

f.  High  intermediate  without  VI 

g.  Low  intermediate  without  VI 

h.  Minimum  without  VI 

2.  Center  of  Gravity 

a.  Forward 

b.  Aft 

3.  Brake  Application  Speed 

a.  +5  knots 

b.  +10  knots 

c.  +20  knots 

d.  + 30  knots 

e.  - 5 knots 
t.  - 10  knots 

4.  Aerodynamics 

a.  No  spoilers 

b.  80%  effective  spoilers 

c.  60%  effective  spoilers 

d.  40%  effective  spoilers 

e.  20%  effective  spoilers 

f.  120%  engine  idle  thrust 

g.  110%  engine  idle  thrust 

h.  90%  engine  idle  thrust 

i.  80%  engine  idle  thrust 

5.  Pilot  Technique 

a.  75%  of  full  metered  pressure 

b.  50%  of  full  metered  pressure 
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Table  4.— Test  Conditions  ( Concluded ) 


Parametric  Studies  (Continued) 

Runway  and  Environmental  System 

1.  Wind 

a.  5 knots 

b.  10  knots 

c.  15  knots 

d.  20  knots 

e.  -5  knots 

f.  -10  knots 

2.  Air  Density 

a.  Hot  day  (83°F/28°C).  hiah  altitude  (5000  ft) 

b.  Cold  day  (-60DF/-51°C) , sea  level 


Landing  Gear  Systems 
1.  Mu-Slip  Curves 

a.  Flat  n~d  peak 

b.  Low  tire  heating 

c.  Tire  inflation  pressure  80?  of  nominal 

d.  Tire  inflation  pressure  120?  of  nominal 


Table  5.— Parameter  Values 


Test  Condition  and 
Parameter  Changed 

Airplane  Model 

B-52 

| KC-135 

F-111 

| Airplane  Parameters 

la  Maximum  Landing  Weight  with  VI 

WA 

300000 

80000 

VI 

267.3 

261 .5 

IYY  x 10  b 

4.927 

.225 

lb  High  Landing  Weight  with  VI 

♦ 

WA 

242500 

68500 

VI  fi 

240.8 

241.5 

IYY  x 10  0 

4.344 

.217 

1c  Low  Landing  Weight  with  VI 

* 

WA 

142500 

52900 

V1  < 

183.3 

211.0 

IYY  x 10‘b 

3.330 

.205 

Id  Minimum  Landing  Weight  with  VI 

★ 

WA 

125000 

48800 

VI  6 

171.1 

202.3 

IYY  x 10  b 

3.153 

.202 

le  Maximum  Landing  Weight  without  VI 

WA  , 

450000 

300000 

80000 

IYY  x 10  6 

8.16 

4.927 

.225 

If  High  Landing  Weight  without  VI 

WA 

370000 

242500 

68500 

IYY  x 10  6 

6.87 

4.344 

.217 

lg  Low  Landing  Weight  without  VI 

WA  , 

245000 

142500 

52900 

IYY  x 10  b 

4.87 

3.330 

.205 

Ih  Minimum  Landing  Weight  without  VI 

WA  * 

20000 

125000 

48800 

IYY  x 10"® 

4.14 

3.153 

.202 

2a  Forward  Center  of  Gravity 

* MAC 

17 

18 

23 

LA 

24.33 

40.89 

19.41 

LB  fi 

25.42 

4.77 

4.98 

IYY  x 10_t> 

4.505 

3.114 

.2035 

* TEST  NOT  COMPATIBLE  WITH  STANDARD  B-52  OPERATION 
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Table  5— Parameter  Values  ( Continued ) 


Test  Condition  and 

Airplane  Model  ! 

Parameter  Changed 

HBPB 

WSBSM 

hhhh 

| Airplane  Parameters 

2b 

Aft  Center  of  Gravity 

% MAC 

35 

35 

60 

LA 

28.46 

44.32 

22.76 

LB  fi 

21.29 

1.35 

1.64 

IYY  x 10'° 

7.515 

4.635 

.221 

3a 

Brake  Application  Speed  + 
5 Knots 

VI 

160.4 

217.4 

228.2 

3b 

Brake  Application  Speed  + 
10  Knots 

VI 

168.9 

225.9 

236.7 

3c 

Brake  Application  Speed  + 
20  Knots 

VI 

185.8 

242.8 

253.6 

3d 

Brake  Application  Speed  + 
30  Knots 

VI 

202.7 

259.7 

270.5 

3e 

Brake  Application  Speed  - 
5 Knots 

VI 

143.6 

200.6 

211.4 

3f 

Brake  Application  Speed  - 
10  Knots 

VI 

135.2 

192.1 

202.9 

4a 

No  Spoilers 

CL 

1.0 

.860 

1 .05 

CD 

.257 

.131 

.144 

4b 

80S  Effective  Spoilers 

CL 

.444 

.4?) 

.434 

CD 

.3082 

.1878 

.2184 

4c 

60X  Effective  Spoilers 

CL 

.583 

.530 

.588 

CD 

.2954 

.1736 

.1998 

4d 

40X  Effective  Spoilers 

CL 

.722 

.640 

.742 

CD 

.2826 

.1594 

.1812 
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Table  5— Parameter  Values  ( Continued ) 


Test  Condition  and 

I” 

1 

Airplane  Model  f 

Parameter  Changeu 

| B-52 

| KC-135 

1 -Ldll 1 

Airplane  Parameters  j 

4e 

20%  Effective  Spoilers 
CL 
CD 

.861 

.2638 

.750 

.1452 

.896 

.1626 

4f 

120%  Engine  Idle  Thrust 
FEO 
KE 

5760 

-8.52 

2880 

-4.26 

108b 
-’  .006 

4g 

110%  Engine  Idle  Thrust 
FEO 
KE 

5280 

-9.59 

2640 

-4.793 

994 
-1  .13? 

4h 

90%  Engine  Idle  Thrust 
FEO 
KE 

4320 

-11.72 

2160 

-5.858 

814 

-1  .384 

4i 

80%  Engine  Idle  Thrust 
FEO 
KE 

3840 

-12.78 

1920 

.6.39 

723 

-1.51 

5a 

75%  of  Full  Metered  Pressure 

- 

- 

- 

5b 

50%  of  Full  Metered  Pressure 

Runway  and  Environmental  System 

la 

5 Knots  Wind 
VW 

8.4 

8.4 

8.4 

lb 

10  Knot  Wind 
VW 

16.8 

16.8 

16.8 

lc 

15  Knot  Wind 
VW 

25.2 

25.2 

25.2 

id 

20  Knot  Wind 
VW 

33.7 

33.7 

33.7 
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Test  Condition  and 
Parameter  Changed 


B-a2 


Airplane  Model 
KC-135  I F-lll 


Runway  and  Environmental  System 


le  - 5 Knot  Wind 
VW 

- 8.4 

- 8.4 

- 8.4 

If  - 10  Knot  Wind 
VW 

- 16.3 

- 16.8 

- 16.8 

2a  Hot  Day 
RHO 

.00189 

.00189 

.00189 

2b  Cold  Day 

.00309 

.00309 

.00309 

RHO 

Landing  Gear  Systems 


la  Flat  fi-a  Peak 

Change  Mu-Slip  Curve 

lb  Low  Tire  Heating 

Change  Mu-Slip  Curve 

lc  Tire  Inflation  Pressure 

Change  Mu-Slip  Curve 

120*  of  Nominal 

Id  Tire  Inflation  Pressure 

Change  Mu-Slip  Curve 

80*  of  Nominal 

2.  PERFORMANCE  STUDIES 


The  performance  studies  provide  a measure  of  the  performance  capability  of  the  brake 
system.  The  tests  performed  fall  into  two  categories.  The  first  defines  the  operation  of  the 
airplane  under  stable  landing  conditions.  The  second  evaluates  the  ability  of  the  brake  system 
to  adapt  to  the  typical  dynamic  operating  conditions  encountered  during  an  actual  landing. 

3.  HYDRAULIC  SYSTEM  STUDIES 

The  hydraulic  system  studies  measure  the  response  of  the  antiskid  valve,  control  box.  and  the 
actual  brake  hydraulic  system.  Specific  tests  were  designed  to  define  both  the  overall  and 
component  performance  of  the  system.  The  results  provide  an  insight  into  aircraft  braking 
system  performance  and  can  be  used  to  further  improve  some  of  the  systems. 

The  test  conditions  as  listed  in  Table  4 were  changed  from  previous  work  (Reference  1 ) to 
provide  a more  comprehensive  and  significant  test  program.  Parameter  variations  were  gener- 
ally made  in  smaller  increments  and  changed  over  a wider  range.  The  numerous  landing  gear, 
pilot  dependents,  and  hydraulic  system  parametric  changes  which  were  included  in  the  pre- 
vious sensitivity  tests  (Reference  1 ) were  excluded,  because  they  are  not  realistic  and/or 
relevant  operational  test  conditions. 
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SECTION  IV 
TEST  RESULTS 


I.  BASELINE 


A baseline  airplane  was  defined  for  each  of  the  three  aircraft  studied.  The  baseline  configura- 
tion is  meant  to  represent  a typical  aircraft  during  its  landing  phase.  The  numerical  values 
of  the  baseline  parameters  are  listed  in  Table  2.  The  values,  when  used  in  the  brake  control 
simulation,  establish  a unique  distance  versus  ground  friction  relation  for  each  airplane. 


i 


Figure  3 relates  the  braking  distance  and  peak  available  ground  friction  (peak  available  mu) 
for  each  of  the  three  baseline  airplanes.  The  data  is  associated  with  the  braking  segment 
only,  the  distances  shown  represent  only  braked  distances,  approach,  flare,  and  transition 
distances  are  excluded.  Peak  available  mu  (such  as  shown  m Figure  3)  is  a computer  input 
that  defines  the  maximum  value  of  friction  available  between  the  tire  and  ground  during  a 
test  condition.  The  value  of  peak  available  friction  was  held  constant  throughout  the  entire 
braking  run  to  generate  the  distance-friction  data  of  Figure  3. 

An  ideal  brake  system  should  operate  at  the  peak  available  friction  value  during  the  entire- 
stop.  For  an  actual  brake  system,  the  instantaneous  coefficient  of  friction  at  the  tire-runway 
interface  depends  on  the  condition  of  the  runway,  tire  properties,  tire  slippage  and  antiskid 
system  characteristics.  Consequently,  the  brake  system  actually  functions  over  a range  of 
friction  somewhat  lower  than  the  peak  value.  Since  the  actual  antiskid  system  has  been  used 
in  the  brake  control  simulation  the  distances  produced  reflect  the  efficiency  of  the  braking 
system. 

Figure  4 is  a plot  of  baseline  braking  distance  efficiency  as  a function  of  peak  available  friction 
for  each  of  the  three  test  airplanes.  Braking  distance  efficiency  is  the  ratio  (in  percent)  of  per- 
fect braking  distance  to  actual  braking  distance  The  efficiency  curves  were  determined  using 
the  data  of  Figure  3.  A detailed  definition  of  braking  distance  efficiency  and  perfect  braking 
distance  is  given  in  Section  V. 

The  KC-135  and  F-l  1 1 efficiency  curves  are  quite  typical  of  second  generation  (early  1960’s) 
brake  control  systems.  The  B-52  efficiency  curve  is  unique  and  requires  explanation.  The 
B-52  brake  control  system  is  a first  generation  (early  1950’s)  system.  The  system  w-as  designed 
as  a tire  saver  (preventing  tire  blow-out  due  to  wheel  lockup)  with  little  attention  to  braking 
efficiency.  The  efficiency  curve  has  three  distinct  regions;  .05  to  .225  mu,  .225  to  .3  mu, 
and  .3  to  ,6  mu.  Each  of  the  three  regions  is  characterized  by  a unique  form  of  brake  control 
operation.  At  values  of  friction  above  .3  the  B-52  is  torque  limited  during  the  entire  braking 
run  (Figure  5a).  Torque  limiting  occurs  when  the  torque  due  to  friction  force  at  the  tire-run- 
way interface  (ground  torque)  exceeds  the  available  brake  torque.  This  is  typical  of  an  under- 
sized brake. 

In  the  region  of  .225  to  .3  mu  a combination  of  brake  control  activity  and  torque  limiting 
occurs  (Figure  5b).  Skidding  activity  occurs  during  the  first  portion  of  the  braking  segment, 
when  insufficient  ground  friction  and/or  vertical  tire  load  exists  so  as  to  cause  torque  limiting. 
As  lift  decreases,  tire  loading  increases  to  the  point  where  ground  torque  exceeds  brake  torque 
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BRAKING  DISTANCE  EFFICIENCY 
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Figure  5.-B-52  H Braking  Activity  Comparison 


resulting  in  torque  limiting  and  the  elimination  of  skid  activity.  The  duration  of  the  initial 
skid  activity  increases  as  ground  friction  decreases,  and  results  in  decreased  braking  distance 
efficiency.  The  decrease  in  efticiency  is  due  to  the  excessively  long  time  between  skids,  when 
brake  pressure  is  zero. 

From  0.05  to  .225  mu  antiskid  activity  occurs  during  the  entire  braking  segment.  In  this 
region  stopping  efficiency  increases  with  decreasing  ground  friction.  This  phenomenon  is  a 
result  of  the  decreased  level  of  pressure  (brake  torque)  that  is  required  to  cause  a tire  to  skid. 
As  friction  decreases  the  system  is  able  to  make  more  efficient  use  of  available  brake  pressuie. 

2,  SENSITIVITY 

The  sensitivity  study  of  the  three  airplanes  involved  changing  a parameter,  or  group  of  param- 
eters, and  observing  the  effect  on  braking  distance.  1 he  braking  distance  associated  with  a 
change  was  then  compared  to  the  baseline  airplane  distance  at  the  same  value  of  peak  avail- 
able mu.  In  this  manner,  the  effect  of  a specific  change  could  be  analyzed  quantitatively. 

The  braking  data  obtained  from  the  brake  control  simulator  represents  an  absolute  distance. 
To  facilitate  the  analysis  of  a parameter  change,  a normalized  distance  has  been  introduced. 

It  is  termed  “percent  baseline  braking  distance”.  The  use  of  a normalized  distance  allows  the 
three  airplanes  to  be  compared  simultaneously.  Percent  baseline  braking  distance  is  defined 
in  Section  V. 

Each  airplane  uses  its  own  baseline  distances  as  normalizing  factors.  Distances  longer  than 
baseline  are  greater  than  100$  and  shorter  distances  are  reflected  as  less  than  100$. 

The  braking  distance  results  from  the  sensitivity  study  have  been  reduced  to  bar  charts, 

Figure  6,  pages  22  through  26.  The  test  conditions  corresponding  to  each  bar  chart  are 
given  to  the  left  and  below  the  chart.  On  each  graph,  the  baseline  braking  distance  percent- 
age for  the  three  airplanes  has  been  plotted  for  four  values  of  peak  available  mu:  0.6,  0.4, 

0. 2,  and  0.1 , (0.15  for  F-l  1 1 ).  During  initial  testing  of  the  F-l  1 1 it  was  found  that  antiskid 
system  operation  on  low  mu  runways  (less  than  .15)  was  sporadic  (See  Volume  II,  Section  X). 
In  an  attempt  to  produce  consistent  data.  .1 5 mu  was  selected  as  the  minimum  runway 
friction  coefficient  for  testing. 

A brief  analysis  of  the  results  follows.  The  raw  data  along  with  the  associated  performance 
indices  can  be  found  in  Volume  II,  Section  VIII. 

A.  AIRPLANE  FLIGHT  CHARACTERISTICS 

1.  Landing  Weight  with  Initial  Velocity  Variation  (Tests  la  through  Id,  Figure  6,  page  22) 

Landing  weight,  initial  airplane  velocity,  and  mass  pitching  moment  were  changed  in  these 
tests.  The  velocity  was  adjusted  to  reflect  changes  in  stall  speed  resulting  from  a landing 
weight  variation.  The  mass  pitching  moment  was  varied  to  reflect  changes  in  load  distribution. 
In  each  test,  brake  application  occurred  0.75  seconds  after  initiation  of  a computer  run.  These 
tests  were  not  applicable  to  the  B-52  since  it  is  not  normal  operating  procedure  to  vary  brake 
application  velocity. 
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1.  LANDING  HEIGHT  WITH  INITIAL  VELOCITY  VARIATION  (Tests  la  through  1b) 

* TEST  NOT  COMPATIBLE  WITH  STANDARD  B-52H  OPERATING  PROCEDURE 


2.  LANDING  HEIGHT  WITHOUT  INITIAL  VELOCITY  VARIATION  (Tests  le through  ih) 

Figure  6.— Sensitivity  Study  Test  Results 
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Figure  6.— Sensitivity  Study  Test  Results  ( Continued ) 
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Figure  6.— Sensitivity  Study  Test  Results  ( Continued ) 
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\s  shown  in  the  bar  charts,  combined  weight,  velocity  and  pitch  moment  changes  resulted 
in  large  braking  distance  variations.  The  distance  variations  result  from  a change  in  the  kinetic 
energy  dissipated  during  braking. 

2.  Landing  Weight  Without  Initial  Velocity  Variation  (Tests  le  thru  111,  Figure  6,  page  22) 

Tests  similar  to  la  thru  Id  were  run  however  the  initial  airplane  velocity  was  held  at  its  appro- 
priate baseline  value.  These  tests  indicate  that  the  airplane  weight  (and  load  distribution)  has 
a large  effect  on  braking  distance. 

3.  Center  of  Gravity  (Tests  2a  and  2b,  Figure  6,  page  23) 

Results  indicate  that  the  aft  location  of  center  of  gravity  causes  a decrease  in  braking  distance. 
This  is  due  to  the  increased  load  placed  on  the  main  gear,  which  increased  the  available  braking 
force.  The  F-l  1 1 appears  to  be  more  sensitive  to  C'G  changes.  This,  however,  is  due  in  part 
to  the  larger  range  over  which  the  F-l  1 1 CG  can  move.  The  B-52  has  braked  nose  and  main 
gears,  consequently,  fore  and  afi  load  transfer  has  no  effect  on  total  available  braking  force 
(as  it  does  with  a brake  main  gear  airplane) 

4.  Brake  Application  Speed  (Tests  3a  thru  3f.  Figure  6,  page  23) 

Brake  application  speed  was  varied  during  these  tests.  Results  indicate  that  additional 
braking  distance  was  required  to  dissipate  the  increased  kinetic  energy  of  the  aircraft. 

5.  Spoiler  Effectiveness  <4a  thru  4e,  Figure  6,  page  24) 

Lift  and  diag  coefficients  were  varied  between  full  landing  spoilers  and  no  spoilers  conditions. 
Braking  distance  increases  as  the  spoilers  become  less  efficient.  The  increased  braking  distance 
results  from:  ( 1 ) loss  of  braking  force  due  to  increased  lift  and  (2)  loss  of  effective  drag 
force  due  to  lower  drag. 

6.  Engine  Idle  Thrust  (Tests  4f  through  4i,  Figure  6.  page  24) 

Increasing  engine  idle  thrust  increases  the  kinetic  energy  to  be  dissipated  during  braking  The 
results  indicate  that  this  resulted  in  longer  stopping  distances. 

7.  Metered  Pressure  (Tests  5a  and  5b,  Figure  6.  page  25 ) 

The  results  from  the  metered  pressure  tests  do  not  show  a general  stopping  distance  trend. 

The  variation  between  aircraft  are  attributable  to  valve  characteristics,  torque  limiting  char- 
acteristics of  the  brake  and  skid  control  system  adaptation  to  conditions.  The  braking  dis- 
tance of  the  F-l  1 1 was  generally  reduced  as  pressure  was  lowered.  The  dynamic  response  of 
the  hydraulic  system  is  slower  at  lower  pressures  reducing  antiskid  cycling,  increasing  system 
efficiency  and  reducing  braking  distance.  The  KC-135.  however,  was  largely  unaffected  by 
metered  pressure  changes  except  at  .6  mu  where  torque  limiting  occurred. 
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Metered  pressure  testing  of  the  B-52  was  not  done  due  to  an  inability  to  obtain  constant 
lower  metered  pressures. 

B.  ENVIRONMENTAL  PARAMETERS 

1. :  Wind  (Tests  la  through  If,  Figure  6,  page  25) 

During  these  tests  the  relative  wind-airplane  velocity  was  maintained,  while  the  airplane-ground 
velocity  and  lift  and  drag  forces  were  changed.  Headwinds  decrease  the  airplane  kinetic  energy 
to  be  dissipated  during  braking  resulting  in  shorter  braking  distances. 

2.  Air  Density  (Tests  2a  and  2b,  Figure  6,  page  26) 

Air  density  was  varied  to  change  the  lift  and  drag  forces  on  the  airplane.  On  hot  days,:  both 
lift  and  drag  decrease,  on  cold  days,  they  increase.  Results  indicate  that  drag  has  the  greater 
effect  on  braking  distance.  On  hot  days,  the  distances  are  longer  because  of  reduced  aero- 
dynamic drag. 

C.  LANDING  GEAR  SYSTEM 

1.  Mu-Slip  Curves  (Test  la  through  Id.  Figure  6,  page  26) 

Variations  in  the  shape  of  the  Mu-slip  curves  were  made  to  analyze  how  the  wheel,  brake  and 
antiskid  system  reacts  to  various  mu-slip  characteristics.  The  randomness  between  the  differ- 
ent aircraft  results  from  the  basic  control  characteristics  of  each  antiskid  system,  however  it 
may  be  generalized  that  braking  distance  decreased  as  the  time  spent  on  the  backside  (negative 
slope)  of  the  mu-slip  curve  increased. 


SECTION  V 

PARAMETER  EVALUATION  CRITERIA 


1.  PERFORMANCE  INDICES 

To  evaluate  the  performance  of  a system,  four  parameters  were  used.  These  performance 
parameters  were 

• Airplane  braking  distance 

• Perfect  braking  distance 

• Braking  distance  efficiency 

• ' "<  baseline  braking  distance 

A.  AIRPLANE  BRAKING  DISTANCE 

Airplane  braking  distance,  X^.  as  measured  on  the  analog  computer  is  the  distance  the  air- 
plane travels  from  brake  application  to  a low-velocity  turn-off  speed  with  the  brake  control 
system  operating. 

B.  PERFECT  BRAKING  DISTANCE 

Perfect  braking  distance,  Xp,  is  the  distance  required  to  stop  an  airplane  when  the  friction 
coefficient  at  the  tire  runway  interface  is  at  its  peak  available  value  during  the  entire  braking 
run.  This  is  the  minimum  distance  in  which  the  airplane  can  be  stopped  at  a given  friction 
value.  Peak  available  friction  is  an  input  to  the  computer  which  s.an  be  a fixed  value  ( .05  to 
.6  typically)  or  variable  (a  function  of  velocity  simulating  a wet  runway). 

Friction  has  been  found  to  be  a function  of  tire  slippage  at  the  tire-runway  intei  face 
(see  Figure  7).  Maximum  braking  is  realized  only  when  the  percent  oi  available  friction  is 
100,  this  occurs  at  about  10%  >hp.  A typical  antiskid  system  allows  slip  to  vary  from  0 to 
100%,  resulting  m fnction  values  les*  than  peak.  Perfect  braking  distance  is  produced  on 
the  computer  by  artificaliy  maintaining  the  percentage  of  maximum  friction  at  100. 

C.  BRAKING  DISTANCE  EFFICIENCY 

Braking  distance  efficiency,  ns,  is  the  ratio  of  ihe  perfect  braking  distance  to  tlm  braked 
airplane  distance. 


ns  * Xp/XA  x 100% 


where:. 


ns  - hrakmg  distance  efficiency 


29 


fcf W«w  w.m ,. 


Xp  = perfect  braking  distance 
= airplane  braking  distance 

Braking  distance  efficiency  indicates  the  degree  ie  which  the  system  meets  its  primary 
requirement  of  stopping  the  aircraft. 

D.  PERCENT  BASELINE  BP  AXING  DISTANCE 

Percent  baseline  braking  distance  is  a dimensionless  distance  ratio  which  facilitates  the 
comparison  of  sensitivity  test  results.  It  is  the  ratio  of  the  braking  distance  resulting  from  a 
parametric  change  to  the  baseline  braking  distance.  The  baseline  braking  distance  of  each 
aircraft  as  a function  of  peak  available  friction  is  given  in  Figure  3. 

2,  SYSTEM  STABILITY 

In  addition  to  the  four  performance  indices,  a measure  was  made  of  the  ability  of  an  antiskid 
system  to  contribute  to  system  stability.  The  criterion  used  to  determine  system  stability 
was  that  the  system  is  termed  unstable  if  the  main  gear  strut  oscillations  diverge  and/or  strut 
oscillations  cause  a reduction  in  brake  pressure.  The  stability  of  the  system  was  measured  by 
determining  the  damping  ratio  necessary  to  cause  instability. 

3.  PARAMETER  RATING  SYSTEM 

The  final  step  in  the  sensitivity  analysis  of  aircraft  braking  performance  was  to  rate  each  param- 
eter. To  facilitate  the  rating  of  a parameter  change  a normalized  distance  referred  to  as  “base- 
line braking  distance  percentage”  was  used.  The  use  of  this  term  allows  the  three  airplanes  to 
be  analyzed  as  a group.  In  order  to  uniformly  and  quantitatively  rate  a parameter  change,  the 
following  formula  was  used: 


n 

2 Percent  baseline  braking  distancej  - 1 00',? 


PR1  = parameter  rating  index 
n = total  number  of  data  points  for  a 
particular  parameter  change 

Percent  baseline  braking  disttMcej  = the 
baseline  braking  distance  percentage 
value  for  the  ith  data  point 


The  parameter  rating  index  (PRl)  as  calculated  above  is  the  average  percentage  deviation  from 
the  baseline  braking  distance.  Thus,  the  vaiue  of  the  PRI  increases  when  a parameter  change 
causes  larger  deviations  from  the  baseline  braking  distance.  The  data  used  in  the  calculation 
of  the  PR!  are  given  in  Volume  II.  Also  included  are  the  final  PRl  values  for  the  dry-stabilized 
landing  conditions. 


SECTION  VI 

SIMULATOR  TO  AIRPLANE  CORRELATION 


Meaningful  results  from  a simulator  can  be  obtained  only  when  the  dynamic  performance 
duplicates  that  of  the  actual  system.  The  brake  control  simulator  used  during  this  study 
is  the  culmination  of  considerable  effort  including  model  development,  subsystem  testing  and 
evaluation,  correlation  with  flight  tests,  and  operational  usage. 

Verification  of  a simulated  aircraft’s  performance  consists  of  comparing  computer  results 
with  results  obtained  from  similar  tests  conducted  during  actual  airplane  flight  tests.  Typical 
data  required  to  simulate  a specific  test  condition  includes  airplane  weight.  C.  G,  location, 
brake  application  speed,  environmental  conditions,  etc.  Key  parameters  within  the  simu- 
lation can  be  adjusted  until  the  desired  level  of  correspondence  is  obtained.  The  following 
items  are  key  points  that  were  considered  during  this  study  in  evaluating  and  obtaining 
simulator-to-airplane  correlation. 

• Stopping  distance 

• Skidding  pressure 

• Number  of  skids 

• Depth  of  skids 

• Rate  in  and  out  of  skids 

• General  control 

No  attempt  was  made  to  duplicate  airplane  stuping  distance  exactly.  Instead,  emphasis 
was  placed  on  producing  the  same  basic  control  characteristics. 

The  airciaft  data  and  flight  test  data  used  for  formulation  and  verification  of  each  airplane 
simulation  was  supplied  by  the  Air  Force.  The  data  provided  for  use  in  the  simulation  was 
adequate,  the  flight  test  data  and  associated  records  used  for  correlation  purposes  were, 
however,  limited.  Ideally,  a number  of  actual  flight  test  records  with  documentation  of 
aircraft  configuration,  initial  conditions,  stopping  distance,  etc.,  are  needed  to  ensure  ade- 
quate airplane  to  simulator  correlation..  In  the  case  of  the  B-52  no  flight  test  records  were 
available,  however,  limited  records  were  provided  on  the  KC-135  and  F-l  1 1 . The  specific 
test  conditions  (KC-135  and  F-l  1 1 ) used  for  correlation  are  listed  in  Table  o.  The  following 
paragraphs  briefly  describe  correlation  procedures  for  each  airplane  simulation. 

1.  B-52 

The  B-52  brake  control  system  posed  unique  problems  during  simulator  setup  and  correlation. 
The  major  concern  during  setup  was  the  locked  wheel  and  skid  detector.  The  detector  is  the 
nuyor  control  element  in  the  brake  system.  The  mechanical  nature  of  the  device  required 
that  it  be  simulated  on  an  analog  computer.  To  correctly  simulate  the  unit,  an  actual  aircraft 
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Table  6.— Correlation  Tests 


Airplane 

Test  program  or 
reference  document 

Test  number 

Test  conditions 

F-111A 

F-111A  Initial  Category  II 
Landing  Performance 
Evaluation 

(See  reference  21) 

Flight  #4 

(Dry 

runway) 

WA  - 53800  VI  = 226.7  fps 

LA  = 20.95  VSTOP  = 38.2  fps 

LB  = 3.44  RHO  = .00222 

FEO  = 950  KE  = 1.28 

Flight  #6 

(Wet 

runway) 

WA  = 55600  VI  = 228 

VSTOP  = 35.3 
LA,  LB.  FEO,  KE,  RHO  UN 
CHANGED  FROM  FLIGHT  #4 

KC-135 

Evaluation  of  a 5-rotor 
Brake  and  Modulated 
Antiskid  System  installed 
on  a KC-135A 
(See  reference  21) 

WA  = 135,200  %MAC  = 19.9 

VI  = 199.62  fps  CL  = .30 

FEO  = 2304  CD  = 155 

KE  = 5.24 

B-52H  flight  manual: 
(See  reference  8) 

T.O  1B-52H- 
1-1 

Figure  A8.1 1 

Case  #1  WA  = 200,000 

Case  #2  WA  = 290,000 

Case  #3  WA  = 370,000 

locked  wheel  and  skid  detector  was  tested  in  the  laboratory  to  determine  its  performance 
A detailed  description  of  the  device,  test  results  and  the  analog  simulation  may  be  found 
in  Volume  II,  Section  111.  In  addition,  stopping  procedure  for  the  B-52  is  unique  and  re- 
quired special  consideration.  The  landing  ground  run  consists  of  two  distinct  segments,  a 
pre-braking  section  and  a braking  section  as  shown  in  Figure  8.  Operating  manuals  foi  the 
B-52H  instruct  the  pilot  to  apply  the  brakes  at  90  knots.  Thus,  the  pre-braking  run  extends 
from  touchdown  to  90  knots.  During  this  segment  aerodynamic  drag  and  rolling  friction 
decelerate  the  airplane  until  the  drag  chute  is  deployed  at  135  knots. 

Flight  test  traces,  pertaining  to  performance  of  the  B-52  brake  control  system,  were  not 
available.  As  a result,  the  basic  control  characteristics  could  not  be  accurately  checked. 
Verification  of  the  simulator  was,  however,  obtained  by  comparison  of  landi  lg  distances 
found  in  the  B-52H  flight  manual  (Reference  8)  with  uistances  generated  by  the  simulation. 
Three  aircraft  weights  were  considered,  a .02  rolling  mu  was  assumed  for  the  prebrake  land- 
ing segment.  Table  7 shows  the  simulator  results  and  the  flight  manual  data.  The  results 
are  not  directly  comparable  (since  the  flight  manual’s  RCR  classification  is  a quahtjt.ve 
measure  of  friction),  however  the  range  and  lower  limit  of  stopping  distance  are  comparable. 

The  general  control  performance  of  the  B-52  antiskid  system  is  typical  of  a first  generation 
skid  control  system.  The  control  performance  is  characterized  by  a series  of  skids  at  a rate 
of  about  1 to  2 skids  per  second.  When  a skidding  wheel  is  sensed,  brake  pressure  is  removed, 
allowing  the  wheel  to  spin-up.  However,  the  wheel  goes  into  another  skid  as  soon  as  brake 
pressure  is  reapplied.  This  skid  cycling  and  prescire  dump-fill  pattern  is  typical  of  the  skid- 
control  operation.  Figures  9 and  10  are  sample  simulator  records  of  antiskid  activity  occur- 
ring during  the  braking  run  on  intermediate  and  low  friction  runways.  It  was  found  that  the 
B-52  brake  becomes  torque  limited  (see  ASD-TR-77-b.  Volume  M.  Section  X)  at  higher  values 
of  friction  coefficient  and/or  aircraft  weight.  This  characteristic  is  i.^ted  in  the  B-52  Flight 
Manual  (Reference  6,  page  7-39).  Torque  limiting  results  in  a constant  stopping  distance  as 
runway  friction  increases.  This  point  is  alluded  to  in  the  B-52  performance  Uu*a  (reference  8), 
when  stopping  distance  becomes  a constant,  even  though  the  RCR  increases. 

2.  KC-135 

Airplane  to  simulator  correlation  was  aided  by  use  of  an  Air  Force  brake  and  antiskid  system 
report  (Reference  15).  This  document  provides  actual  flight  test  data  and  sample  records  of 
wheel  and  brake  pressure  time  histories  for  both  wet  and  dry  runways.  !t  must  be  noted  that 
numerous  inconsistencies  exist  in  the  time  histories.  Specifically  the  wheel  speed,  antiskid 
voltage  and  brake  pressure  for  the  wet  runway  conditions  do  not  appear  to  be  compatible 
with  one  another.  Additionally,  the  antiskid  voltage  level  m both  the  wet  and  dry  cases  is 
not  consistent  with  previous  experience  and  current  results.  Consequently,  only  the  brake 
pressure  and  wheel  speed  time  histories  along  with  the  braking  distances  reported  in  Ref- 
erence 15,  and  previous  experience  with  the  Mark  il  antiskid  system  (found  on  the 
KC-135)  were  used  to  evaluate  simulator  correlation. 

A portion  of  flight  number  2-7,  dry  runway,  ( Reference  1 5)  is  shown  in  Figure  1 1 . along 
with  an  equivalent  condition  obtained  from  the  simulator.  The  skid  control  performance 
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Figure  8.—B-52H  Stopping  Procudure 


■Simulator  to  Airplane  Correlation  Results,  B-52 
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STOPPING  DISTANCE  FROM  TOUCHDOWN 
RUNWAY  CONDITION  READING 

DATA  COMPILED  FROM  TO  1B-52H-1-1  FIGURE  A8-11 

COMPUTER  GENERATED  DISTANCE  ASSUMED  A ROLLING  MU  OF  .02  PRIOR  TO  BRAKING. DURING  BRAKING 
MU  IS  AS  SPECIFIED  IN  THE  TABLE 
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Figure  9.—B-52H  Simulator  Performance  (Baseline  Airplane) 
Intermediate  Mu  (.25) 


is  characterized  by  skidding  activity  at  regular  intervals.  The  skid  rate  is  about  one  per 
second.  Between  skid  cycles  pressure  is  ramped  on  by  the  system  at  a rate  of  about  200  to 
250  psi/sec.  Skidding  occurs  at  approximately  600  psi.  Both  the  simulator  and  flight  test 
data  show  these  characteristics.  The  aii plane  braking  distance  was  1750  feet,  while  the 
distance  obtained  on  the  simulator  was  1772. 

3.  F-l 1 1 


The  F-l ! 1 simulator  performance  was  compared  with  flight  tests  4 through  7 (Reference  21 ). 
Portions  of  tests  4 and  6 are  shown  in  Figures  12  through  14,  along  with  an  equivalent  con- 
dition obtained  from  the  simulator.  It  should  be  noted  that  the  wheel  speed  signals  in  the 
flight  test  traces  do  not  reflect  skidding  activity  consistent  with  the  antiskid  signal  and  brake 
pressure  traces.  The  wheel  speed  traces  would  appear  to  be  heavily  filtered.  It  is  a common 
instrumentation  practice  to  filter  signals  to  reduce  noise.  However,  this  practice  also  atten- 
uates the  amplitude  of  a rapidly  changing  signal  (such  as  a skid).  A comparison  of  the  traces 
in  Figures  I 2 and  13  shows  that  the  skid  and  recovery  pressure  levels,  antiskid  value  voltage 
levels  and  general  control  characteristics  of  the  actual  aircraft  are  reproduced  by  the  simulator. 
In  the  dry  runway  example,  major  skid  cycling  occurs  at  a rate  of  about  one  skid  every  1 to 
2 seconds;  in  addition,  a higher  frequency  skid  cycling  (about  10  cps)  has  also  been  repro- 
duced, Figure  13.  Resulting  skid  cycling  pressures  are  about  1200  to  1400  psi.  The  actual 
stopping  distance  of  flight  4 was  1 753,  while  the  corresponding  simulator  distance  was  1 748. 
Pressure  and  antiskid  signal  voltage  levels  and  cycling  rate,  are  reproduced  in  the  wet  runway 
condition  (flight  6)  shown  in  Figure  14.  The  simulator  braking  distance  for  this  condition 
was  6006  feet,  while  the  aircraft  distance  was  5956. 
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SECTION  VII 

DETERMINATION  OF  SIGNIFICANT  PARAMETERS 


1.  PARAMETER  RATINGS 

The  parameter  rating  index  (PR1)  was  used  to  rank  the  parameter  changes  according  to  their 
effect  on  braking  distance.  The  value  of  the  PRI  is  the  average  percentage  deviation  from  the 
baseline  braking  distance.  Thus,  the  larrsr  the  PRI.  the  greater  impact  the  parameter  lias  on 
•airplane  braking  distance.  Based  on  the  PRI,  the  parameter  changes  have  been  arranged  in 
numerical  order,  and  the  results  are  listed  in  Table  8.  The  table  presents  the  average  PRI  for 
the  three  aircraft  tested.  Although  the  tire-ground  friction  coefficient  (mu)  is  a predominant 
influence  on  airplane  braking  performance,  it  is  not  included  in  these  tables  because  the  PRI 
rating  methodology  was  used  to  determine  variable  significance  for  a range  of  mu  values. 

2.  SIGNIFICANT  PARAMETERS 

The  second  step  in  the  rating  of  parameters  was  to  determine  which  parameter  changes  1 ave 
a significant  effect  on  stopping  distance.  It  was  decided  to  consider  all  parameters  having  a 
PRI  greater  than  2.0  as  being  significant..  A value  of  2.0  represents  a 2%  change  in  the  base- 
line braking  distance.  The  repeatability  fo  the  analog-hardware  simulation  itself  results  in 
1%  variations. 

Based  on  the  above  criteria,  Table  8 has  been  reduced  and  peak  available  mu  included.  In 
addition,  the  parameter  changes  have  been  summarized  by  combining  related  tests  under  a 
general  heading.  The  resulting  list  of  parameters  having  a significant  effect  on  braking  dis- 
tance is  given  in  Table  9. 
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Table  8—  Parameter  Rating 
Composite  F t 11,  B-52,  KC-135 


RANK 

TEST 

CONDITION 

DESCRIPTION 

mmm 

PR  I 

1 

la 

MAXIMUM  LANDING  WEIGHT  WITH  VI 

43.88 

2 

Id 

20  KNOT  WIND 

37.67 

3 

3d 

BRAKE  APPLICATION  SPEED  + 30  KNOTS 

35.92 

4 

4a 

NO  SPOILERS 

33.50 

5 

1c 

15  KNOT  WIND 

29. 75 

6 

le 

MAXIMUM  LANDING  WEIGHT  WITHOUT  VI 

23.42 

7 

If 

-10  KNOT  WIND 

23.33 

8 

3c 

BRAKE  APPLICATION  SPEED  + 20  KNOT 

23.25 

9 

4e 

20*  EFFECTIVE  SPOILERS 

22.42 

10 

lb 

10  KNOT  WIND 

20.50 

n 

5b 

50*  OF  FULL  METERED  PRESSURE 

15.71 

12 

lb 

HIGH  INTERMEDIATE  LANDING  WEIGHT  WITH  VI 

15.25 

13 

4d 

40*  EFFECTIVE  SPOILERS 

14.95 

14 

Id 

MINIMUM  LANDING  WEIGHT  WITH  VI 

13.50 

lh 

MINIMUM  LANDING  WEIGHT  WITHOUT  VI 

13.50 

16 

If 

HIGH  INTERMEDIATE  LANDING  WEIGHT  WITHOUT  VI 

13.08 

17 

3f 

BRAKE  APPLICATION  SPEED  -10  KNOTS 

12.75 

18 

3d 

BRAKE  APn’  ’CATION  SPEED  + 10  KNOTS 

11.83 

19 

le 

-5  KNOT  •» 

11.33 

20 

la 

5 KNOT  WIND 

11.17 

21 

lc 

LOW  INTERMEDIATE  LANDING  WEIGHT  WITH  VI 

10.63 

22 

4c 

60*  EFFECTIVE  SPOILERS 

8.83 

23 

la 

FLAT  MU-SLIP  PEAK 

8.58 

24 

2b 

COLD  DAY 

8.17 

25 

5a 

75*  OF  FULL  METERED  PRESSURE 

8.14 

26 

3e 

BRAKE  APPLICATION  SPEED  - 5 KNOTS 

7.00 

27 

LOW  INTERMEDIATE  LANDING  WEIGHT  WITHOUT  VI 

6.75 

28 

lb 

LOW  TIRE  HEATING 

6.58 

29 

2a 

FORWARD  CENTER  OF  GRAVITY 

6.50 

30 

lc 

TIRE  INFLATION  PRESSURE  80S  OF  NORMAL 

6.33 

31 

3a 

BRAKE  APPLICATION  SPEED  + 5 KNOTS 

5.82 

32 

2a 

HOT  DAY 

5.75 

33 

Id 

TIRE  INFLATION  PRESSURE  120S  OF  NORMAL 

5.42 

34 

4i 

80S  ENGINE  IDLE  THRUST 

5.25 

35 

4b 

80S  EFFECTIVE  SPOILERS 

3.75 

36 

2b 

AFT  CENTER  OF  GRAVITY 

3.67 

37 

4f 

120*  ENGINE  IDLE  THRUST 

3.58 

38 

4h 

90S  ENGINE  IDLE  THRUST 

2.75 

39 

4g 

110S  ENGINE  IDLE  THRUST 

2.00 
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Table  9.— Significant  Parameters 


Parameters  affecting  braking  distance  by  more  than  2% 
(Listed  in  decreasing  order  of  significance) 

1.  Peak  available  ground  friction 

2.  Landing  weight  with  initial  velocity  changes 

3.  Wind 

4.  Brake  application  speed 

5.  Spoiler  effectiveness 

6.  Landing  weight  without  initial  velocity  changes 

7.  Metered  pressure  level 

8.  Mu-slip  curve  shape 

9.  Air  density 

10.  Fore  and  aft  center  of  gravity 
11..  Engine  idle  thrust 
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SECTION  VIII 

SELECTION  OF  PERTINENT  PARAMETERS 


Figure  1 5 is  a How  chart  where  each  block  represents  a major  step  of  analysis  in  the  formula- 
tion of  the  prediction  equation.  Figure  16  shows  the  breakdown  of  the  entire  Task  II. 

Table  10  lists  those  significant  parameters  remaining  after  excluding,  from  Table  9,  the  param- 
eters that  depend  on  pilot  technique,  or  are  design  constraints  and/or  are  outside  the  scope  of 
present  work.  The  following  paragraphs  present  the  reasoning  for  further  refinement  to  the 
list. 


1,  PEAK  AVAILABLE  MU  AND  MU-SLIP  CURVES 

Parameter  number  1 1 in  Table  9,  namely  mu-slip  curve  shapes,  are  not  independent  variables 
and  therefore  are  inherent  parts  of  or  dependent  upon  the  peak  available  mu  at  the  tire-runway 
interface.  One  of  the  a priori  requirements  when  forming  dimensionless  groups  is  that  each 
term  or  group  be  independent.  Therefore,  Mu-slip  curve  shapes  are  not  listed  on  Table  10. 


2.  AERODYNAMIC  LIFT  AND  DRAG 

Lift  and  drag  coefficients  are  interdependent  variables.  Therefore,  the  ratio  C'l/Cq  was 
chosen  as  an  independent  variable.  Wherever  applicable,  drag-chute  caused  drag  was  added 
to  the  aerodynamic  drag  and  the  sum  treated  as  the  total  aerodynamic  drag. 

3,  HEAD  OR  TAIL  WIND 

In  the  previous  sensitivity  study  report,  reference  (1),  it  was  pointed  out  that  a single  curve 
would  suffice  to  describe  the  variation  of  both  wind  velocity  and  brake  application  speed. 
Thus,  a five  knot  wind  velocity  was  the  same  as  a five  knot  change  in  the  brake  application 
velocity.  This  comparison  was  made  for  the  aircraft  considered  in  present  study  and  the 
results  are  shown  in  figures  17,  18,  and  19,  Clearly  the  aforementioned  axiom  does  not 
apply  at  all  for  the  B-52  and  holds  only  at  high  mu  (0.4/i  and  higher)  conditions  for  the 
KC-135  and  F-l  1 1 airplanes.  Although  the  above  axiom  does  not  apply  in  this  case,  the 
wind  component  can  be  converted  to  a velocity  change  and  therefore  cannot  be  considered 
us  an  independent  variable  apart  from  the  brake  application  velocity,  A correction  factor  is 
therefore  needed  to  convert  a given  wind  component  into  an  equivalent  velocity  component 
to  be  used  in  the  prediction  cquation/model.  These  correction  factors  were  computed,  based 
on  the  sensitivity  study  data  and  are  discussed  in  detail  in  Appendix  A. 

4,  LANDING  WEIGHT  AND  BRAKE  APPLICATION  VELOCITY 

Nominally,  the  touchdown  speed  and,  consequently,  the  brake  application  speed  are  a direct 
function  of  gross  weight;  the  higher  the  weight  the  higher  the  touchdown  speeds.  The  pilot 
technique  and/or  training  (e.g.,  see  the  later  paragraph/discussion  on  B-52  operations)  may, 
however,  result  in  situations  where  the  pilot  applies  brakes  at  the  same  speed  regardless  of  the 
gross  weight  or  touchdown  speed.  It  is  thus  interesting  to  study  the  effect  of  weight  variation 
(without  velocity  variation)  upon  the  stopping  distance. 


Figure  15.— Block  Diagram  for  Task  II  Analysis 


Figure  1 6— Sensitivity  Analysis  and  Braking  Prediction  Subsystem  (BPSS) 


Table  10.— Significant  Parameters 


1. 

Peak  available  mu 

2. 

Drag  device  effectiveness 

3. 

Head  or  tail  wind 

4. 

Brake  application  speed 

5. 

Landing  weight 

6. 

Air  density 

7. 

Engine  idle  thrust 

8, 

Center -of-gravity  location 
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Figure  17.— Effect  of  Wind  and  Brake  Application  Velocities  on  Stopping  Distance . B-52 


A USAF  contract  study,  reference  (24),  had  shown  that  for  the  F-4  airplane  this  effect  was 
less  than  one  percent  in  terms  of  stopping  distance,  in  the  previous  sensi'ivity  study,  refer- 
ences ( 1 1 and  (2), .the  tests  on  the  B-747  simulator  showed  this  effect  to  be  less  than  two 
percent.  The  matter  was  therefore  not  pursued  any  further  at  that  time,  as  the  criteria  for 
selecting  significant  parameters  was  to  detect  a change  in  the  stopping  distance  of  more  than 
two  percent. 

In  the  current  study  it  was  noted  that  the  flight  manuals  for  the  B-5 2.  reference  (8),  instruct 
the  pilot  to  apply  brakes  at  the  same  speed  regardless  of  the  landing  weight  or  touchdown 
speed,  e.g.,  90  knots  for  B-52  G/H  models  (70  knots  for  earlier  versions).  This  necessitated 
generation  of  B-52  simulator  data  in  a fashion  so  as  to  duplicate  the  actual  aircraft  operation 
It  was,  however,  decided  to  obtain  simulator  data  for  all  three  airplanes  (B-52,  KC-135  and 
F-l  1 1 ) to  study  the  weight-alone  variation  effect.  These  data  have  been  plotted,  figures  20 
to  22,  as  percent  change  in  weight,  7c  Aw,  versus  percent  change  in  braking  distance, '/  As; 
the  change  being  calculated  as  ± percentage  values  from  the  defined  baseline  (see  section  11 
for  baseline  definition). 

The  data  for  the  B-52,  figure  20,  is  consistent,  except  for  two  data  points  at  0.2  mu.  It 
shows  a direct  relationship  between  weight  and  distance;  an  increasing  weight  resulting  in 
higher  braking  distance  and  vice  versa.  The  deviations  at  0.2  mu  are  due  to  torque  limiting 
operation  as  well  as  the  skid  control  system  transients  at  that  particular  friction  level.  This 
will  become  more  clear  when  mu-braking  efficiency  plots  are  discussed  in  a later  section  (see 
section  X). 

The  data  for  KC-1  35  and  F-l  1 1 , figures  21  and  22,  show  a mixed  trend  in  that  an  increase 
in  weight  does  not  necessarily  cause  an  increase  in  distance  and  vice  versa.  The  skid  control 
system  efficiency,  available  friction  level,  the  landing  gear  foot  print  (geometry),  the  trans- 
fer/distribution of  airplane  weight  among  the  gears  and  the  strut  stability,  could  individually 
or  collectively  influence  the  resulting  braking  distance  from  one  landing  weight  to  another. 
The  data,  however,  can  be  grouped  together  into  shaded  areas  and  certain  trends  established 
when  the  torque  limited  braking  cases  (designated  TL  in  the  figures)  are  excluded. 

It  was  explained  in  reference  (1 ),  page  78,  as  to  why  only  brake  application  velocity  and  not 
both  the  landing  weight  and  velocity  were  chosen  as  significant  independent  parameters.  In 
the  case  of  B-52,  however,  a special  handling  of  the  weight  variation  (without  velocity  vari- 
ation) data  was  necessary  as  weight  does  not  appear  as  a separate  variable  in  the  nondimen- 
sionai  terms.  An  equ:''alent  velocity  change  was,  therefore,  calculated  for  each  weight 
change.  This  will  be  demonstrated  by  an  example  in  Appendix  B.  The  weight-alone  variation 
data  for  KC-135  and  F-l  1 1,  figures  21  and  22,  were  not  included  in  the  prediction  model 
analysis  as  the  data  for  weight  and  velocity  varying  together  (normal  approach)  was  also 
obtained  on  the  simulator  and  accordingly,  used  in  the  analysis.  Due  to  the  unique  B-52  oper- 
ating procedures  of  delaying  brake  application  until  90  knots,  the  rolling  distance  or  the  pre- 
braking distance  from  touchdown  to  brake  application  amounts  to  a considerable  value,  and 
should  be  calculated  separately  and  added  to  the  braking  distance  to  obtain  the  needed  run- 
way length/ground  roll.  Typical  prebraking  roll  distances  were  calculated  for  the  B-52  and 
are  shown  in  figure  23.  The  mathematical  expression  used  for  these  calculations  is  explained 
in  Appendix  C. 


-r  — ***** 


ASSUMPTIONS: 

1 Rolling  V ■ 0.02 

2 Drag  chute  deployment  speed  ■ 135  Knots 


3 CD  . .152  for  VToucMoim  to  135  Knots 

4 CQ  * .321  for  135  Kn  to  90  Kn 


5.  CENTER-OF-GRAVITY 

The  CG  location  effect  considered  in  the  analog  simulation  was  of  the  form: 

p*  LA 

* “ LA  + LB  + p-HB 

where  LA,  LB,  and  HB  are  geometric  distances  (see  Table  2),  and  p is  the  coefficient  of 
available  friction.  The  parameters'!'  and  (l-*) thus  determine  the  respective  loads  carried 
by  the  main  and  nose  gears  for  given  c.g.  locations.  Because  the  CG  variation  of  any  correla- 
tion ground  vehicle  would  be  minimal,  and  because  geometric  similarity  with  an  airplane 
would  be  almost  impossible  to  achieve,  it  was  decided  to  consider  only  the  coefficient  of 
available  friction  as  the  independent  variable. 

6.  OTHER  PARAMETERS 

The  remaining  parameters  (p  and  Fe)  are  independent  variables  and  require  no  discussion. 
From  the  preceding  paragraphs,  it  follows  that  the  pertinent  variables  are. 


Braking  distance 

(s) 

Available  mu 

(p) 

CL/CD  ratio 

(CL^CD) 

Brake  application  speed 

(v) 

Air  density 

(p) 

Engine  idle  thrust 

<Fe) 

where  (s)  is  th  * dependent  variable  and  all  others  are  independent  variables. 


SECTION  IX 

DEVELOPMENT  OF  PREDICTION  MODEL 


1,  COMPONENT  EQUATIONS 

The  first  step  in  forming  a prediction  model  is  to  identify  the  pertinent  and  independent 
variables.  This  step  is  by  far  the  most  important  because  the  validity  of  the  results  depends 
on  the  correctness  with  which  the  pertinent  factors  are  selected.  For  this  study,  as  explained 
in  Section  VIII.  This  required  the  combining  of  some  of  the  interdependent  variables  listed  in 
Table  10.  The  list  of  resultant  independent  variables  is  shown  in  Table  1 1 . 

The  second  step  is  to  express  the  secondary  quantity  (dependent  variable)  as  a function  of 
the  primary  quantities  (independent  variables),  so  that: 

s = F(g,  v,  p,  Fe,  p,  CL/CD)  (1) 


where- 


s = braking  ;top  distance 
g = acceleration  caused  by  gravity. 

The  dimensional  matrix  that  can  be  formed  for  the  fundamental  units  (mass,  length,  and  time) 
of  the  seven  parame'ers  in  Equation  1 is  of  rank  3,  so  that,  according  to  Buckingham’s  it 
theorem,  these  would  yield  four  independent  7T  terms.  By  inspection  and  analysis,  they  can  be 
written  (sg/v2),  (p),  (Cp/Cp),  and  (pv6/Feg2).  Thus: 


(sg/v2)  = F (p,  Cp/Cp,  pv^/F0g2)  (2) 

or.  , ir  j ) = F (7T-t,  tr^,  7r^)  *2a) 

where. 

tr  | = sg/v2 
7i-i  = p 

713  = CL/Cp 
7r4  = pvbl Fe  g2 

Appendix  C,  reference  ( 1 ).  shows  the  detailed  analysis  of  arriving  at  Eq.  2 and  Eq.  2a. 

Section  XII,  Volume  II  of  this  report*  ASD-TR-77-6)  shows  the  calculations  of  7r  terms 
(numerics)  using  raw  data  from  Task  HA  simulation.  The  application  of  dimensional  analysis, 
including  the  pi  theorem,  leads  to  a type  of  equation  involving  an  unknown  function,  of  which 


Table  1 1—  Pertinent  Independent  Variables 


Variable 

Notation 

Available  mu . . 

P 

Brake  application  speed  ..............  ..  ..  ...  . 

V 

Drag  device  effectiveness  

cL/cD 

Engine  Idle  thrust 

Fe 

Air  density  

P 
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Eq.  2 is  an  example.  Before  a prediction  equation  can  be  formulated,  the  nature  of  the 
function  must  be  determined.  This  cannot  be  accomplished  by  dimensional  analysis,  but 
it  can  be  done  from  analysis  of  laboratory  observations. 


The  best  procedure  for  evaluating  a function  is  to  arrange  the  observations  so  that  all  but 
one  of  the  pi  terms  containing  the  independent  variables  in  the  function  remain  constant. 

Then  the  remaining  independently  variable  pi  term  is  varied  to  establish  a relationship 
between  it  and  the  dependent  variable  (tr  j term).  Section  XIII,  Volume  11  of  this  report 
(ASD-TR-77-0)  shows  this  arrangement  of  experimental  observations  (tr  terms)  for  all  three 
airplanes  under  consideration.  This  procedure  is  repeated  for  each  of  the  pi  terms  in  the 
function:  the  resulting  relationships  between  tt]  and  the  other  individual  pi  terms  are  called 
component  equations.  Statistical  curve  fitting  computer  programs  were  used  to  generate  the 
component  equations  (see  Appendix  D).  A summary  of  the  equations  is  listed  in  Table  1 2, 

It  should  be  noted  that  m some  of  the  component  equations,  the  exponent  for  the  (714) 
term  has  been  modified  by  a Ap  term  where  A p equals  Ip(std-day)  - p(non-std-day)]  in 
lb-sec~/ft^,  the  terms  standard  day  and  non-standard  day  are  defined  in  Appendix  E.  Tins 
was  necessary  to  arrive  at  a satisfactory  relationship  between  (7rj)  and  (*4)  where  satisfactory 
implies  an  acceptable  correlation  error  between  the  actual  and  predicted  values,  e.g..  ± 5 7c.  This 
modification  of  the  (^4)  exponent  was  needed  at  all  p conditions  for  the  B-52  and  only  at  0.2p 
or  less  for  the  KC-135  and  the  F- 1 i 1 airplanes.  The  mechanism  for  arriving  at  the  numerical 
value  of  the  modifier  (Ap  term)  is  explained  by  an  example  calculation  in  Appendix  E. 


An  examination  of  the  previous  sensitivity  study,  reference  ( 1 ) showed  that  even  thougli 
air  density  (p)  lad  been  included  as  an  independent  variable  (as  part  of  the  ir 4 term)  the 
density  change  data  collected  on  the  simulator  had  inadvertently  been  left  out  when  formu- 
lating component  and  prediction  equations  in  the  numerical  form.  This  caused  some  concern 
as  to  the  accuracy  of  the  previously  reported  prediction  equations  for  the  Boeing  727,  737, 
747  and  the  USAF  C-141  and  F-4  airplanes,  reference  ( 1 ).  All  the  component  equations  and 
the  corresponding  prediction  equations  were,  therefore,  recalculated  by  including  the  data 
points  for  density  variation.  A check  of  the  correlation  errors  showed  that  a modified  (774) 
exponent  would  be  needed  for  F4  at  all  p conditions  and  only  for  wet  runways  for  the  other 
four  airplanes.  The  accuracy  of  the  remaining  prediction  equations  and  their  application 
range  were  no*  affected.  Accordingly,  the  revised  prediction  equations  with  mollified 
exponents  are  shown  in  Table  13. 

2.  GENERALIZED  FUNCTIONS 


When  the  component  equations  have  been  determined,  they  are  combined  in  a certain  man- 
ner to  give  a general  relationship.  It  is  possible  for  some  of  the  component  equations  to  be 
combined  by  multiplication,  while  others  require  addition  in  die  formation  of  the  resultant 
prediction  equation.  In  general,  these  two  methods  are  adequate  for  the  majority  of  engineer- 
ing problems.  For  the  stopping  distance  problem,  the  analysis  showed  that  the  prediction 
equation  should  be  formed  by  multiplication.  The  necessary  and  sufficient  conditions  to  be 
met  lor  the  Junction  to  be  a product  were  developed  and  translated  into  tests  of  validity.  A!! 
aspects  of  the  development  of  prediction  equations  discussed  in  this  paragraph  are  detailed 
in  Appendix  E,  reference  ( 1 (.  The  major  equations  of  interest  arc  repeated  in  succeeding 
paragraphs. 
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Table  12— Summary  of  Component  Equations 


Airplane 

model 

u* 

Equation 

Eo. 

No. 

B-52 

0.6 

(n,)  = 2.3931  (ir2)  ■•040° 

(3) 

(*,)  - 2.4772  (.j)  -0^-OVm  SP- 

(4) 

<»,)  • 4.4960  (*4)  9-054c  ' •°68°3 

(5) 

0.4 

( tt 1 ) = 2.3931  (tt2)  '-04007 

(3) 

(.,)  * 2.4815  (.3)  -1285  - -2227SSP 

(6) 

4.0297  <*4>  8.484P-.05S4 

(7) 

0.23 

/ \ f-n  \ “2.2504 

(ffj ) - .1 747  (^2) 

(8) 

(.,)  - 4.8081  (.3)  -'3006  - .0774USP 

(9) 

(tr1 ) - 24.5495  (*4)  17*75Ap  ' *1836 

(10) 

0.2 

(1^)  = 2.3637  (tt2)  ‘-5166 

(ID 

(„,)  . 5.4324  (,3)  -090S1-.07014XSP 

(12) 

, \ oc  ,7yl0  . > 12.5Ap  - .1803 

( tt-j ) = 26.3742  (*4) 

(13) 

0.1 

(tt1  ) = 2.3637  (tt2)--5166 

(ID 

(»,)  - 7.9153  (,3!  -1H44-.13.WSP 

(14) 

(,,)  = 43.8256  (»4)  16'"4p  - ',962 

, 

(15) 

‘Baseline  value  of  p used  in  the  data  set 
“%  SP  is  the  percentage  of  the  spoiler  configuration. 


Table  12— Summary  of  Component  Equations  ( Continued ) 


Airplane 

model 

u * 

Equation 

Eq- 

No. 

KC-135 

0.6 

(»,)  = .71355  (*2)  -*9420 

(16) 

(.,)  = 1.1046  (.3)  -19028-.10895JSP 

(17) 

("■,)  = 2.1672  (tt4)  ‘*05414 

(18) 

0.4 

(*,)  = .71355  (t)2)  -‘9420 

(16) 

(,,)  = 1.68155  (,3)  -22326-. 2071 8*SP 

(19) 

(ir, ) = 4.39388  (tt4)  '-08228 

(20) 

i 

0.2 

(*-,)  = .5893  (tt2)  '1'1410 

(21) 

j 

(„,)  . 3.38793  (,3)  -25766  -02665" 

(22) 

(it, ) = 22.44915  (^4)  3,5Ad  '-15208 

(23) 

0.1 

(it,)  = .5893  (tt2)  -1*1410 

(21) 

(„,)  - 8.0898  (,3)  ',485-'6«%SP 

(24) 

(Tr-j)  = 82.26867  (Tr4)  12,65Ap  “ -200863 

(25) 
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Table  12.— Summary  of  Component  Equations  (Concluded) 


u teggmxssaaBmm — r ... 


Airplane 

model 

u* 

Equation 

R 

F-lll 

.6 

(tt]  ) = .7760  (n2)  '1J411 

(26) 

(tt1  ) = 1.3973  (tt3)  *1890  --1286%SP 

(27) 

(tt1  ) = 2.9702  (tt4)  --05843 

(28) 

.4 

(*,)  = .7760  (ir2)  '11411 

(26) 

(ir,)  - 2.1722  (n3)  * 2060- . 1 448%SP 

(29) 

(rr1  ) = 5.1797  (*4)  ‘•°6742 

(30) 

.2 

fy)  = .9467  (tt2)  ’1'0490 

(31) 

(»,).  5.2212  <»3>  -,885-'185ISP 

(32) 

(tt]  ) -•  42.8987  (*4)  *-1648 

(33) 

.1 5 

(ir1 ) = .9467  (tt2)  '1,049° 

(31) 

(w1 ) = 6.6120  (tt3)  *1656- ,C698%SP 

(34) 

(w1 ) = 66.2844  (n4)  8,8Ao  ' ,1797 

-i 

(35) 

i * 

I 

i 

s 

\ ' 

X 

i 

t 

i 

l 

1 

5 


i 

i 

j 

! 

s 

I 

1 


Table  13.— Revised  Prediction  Equations  for  727,  737,  747,  C-141  and  F-4  Airplanes 


Revised  Prediction  Equations  for  F-4 


Equation 


Oi,)  - 1.7679  (»2)-°*9236  (^3)  C.3485-.1185X  spJ  ^ [6.1  4? -.0468] 
(#j ) * 2.2932  (Wj)  Co.5259-.0783X  Sp]  ^ [l2.26  A? -.0649] 

(»,)  ■ 21 .5085  (.j)"'-'580  (,3)  t«»-l>039*  Sp]  ^ [lZ-8U?-.2419] 


* The  original  prediction  equations  are  listed  in  Reference  (1  ),  Table  14,  p.92 


Revised  Prediction  Equations  • Wet  Runways 


Air- 

plane 

H 

Equation 

u 

EM 

727 

.167 

(*1)  « 1.9647  Gf2)-1*16  (ir3)  L-l^-.OOlX  SP]  ^ Bs.l8-4f-.1108] 

(14) 

737 

.141 

(*,)  « 0.6595  (Tr2)"-7626Cw3)  C2344-.1271*  sp](*4)  0 8.36-4?  *.01 60] 

(15) 

747 

.125 

(w1)  * 0.6704  (w2  )"*8104(ir3)  L3568-.07X  SP]  ^ (21.6-4?  -.0143] 

(16) 

C-141 

.2251 

(•j)  - 2.4953  (ir2)'1*0326(ir3)  £<«69-G002X  S^}  03.48  4? -.0979/ 

07) 

F-4 

.278 

(*,)  « 4.6672  (w2rlj665  (w3)  C.7498-.0257X  SP]  Q 8.29  ^f-.l  320 

....  . ..  J 

(18) 

When  the  component  equations  (see  Table  i 2)  are  combined  by  multiplication,  the  predic- 
tion equation  is  of  the  form: 

JT)  =(C)(7Tj)  3.4 (Wj ) 2^ 4 (T  j)  2,3 
where  the  bar  denotes  a constant  (held)  value. 

The  analysis  shows  that  the  value  of  the  constant  term  C is  of  the  form 

C = — — — — - — — 2 (37) 

(F  (ff2>  ff3’  n4>‘ 

Thus  the  prediction  equation  is  of  the  form: 


F (*2-  tr3, 7r4)  = 


F Qr2,  i3,  tt4)  F (tt2.  tt.i,  tt4)  F (ftp,  if 3,  tt4) 
{F  (5r2,  W3,  7t4)J 2 


(38) 


The  equations  constituting  a test  for  the  validity  of  Eq.  38  are  shown  to  be  (see  Appendix  E, 
Reference  ( 1 ). 


F (ir2,  w3,  7t4)  F (ir2,  rr3,  7r4)  F (jt2,  jt3,  ir4)  F (ir2>  ?3,  rr4) 
[F(?2,f3,f4)l2  lF(72,ff3,¥4)]2 

or  F (ir2,  ?4)  F(¥2,¥3,ir4)  F (?2,7r3,i4)  F(?2- ^3.^4) 

IF  (?2,?3,7r4)]2  . (F  (i2,  ff3, 7r4)j 2 


(39) 


(39a) 


The  values  and  ir3  are  values  of  tt2  and  ff3  held  constant  at  some  value  other  than  f2 
and  f3.  Thus  from  the  observed  data:. 


ir2  = 0.6 
f = 0.4 

#1  = 0.2 

etc. 


the  primary  set  of  data,  for  example 
supp’“mentary  sets  of  data 


If  the  supplementary  sets  of  data  satisfy  either  Eq.  39  or  39a,  the  general  equation  can  be 
formed  by  multiplying  the  component  equations  together  and  dividing  by  the  constant,  as 
indicated  in  Eq.  38, 

This  test  was  applied  to  all  available  data  (component  equations);  the  results  are  shown  in 
Table  14  clearly  indicating  the  validity  of  the  approach.  Table  D-l  contains  the  details  of 
this  calculation. 


Another  test  of  validity  was  to  calculate  the  value  of  the  constant  term  C of  Eq.  36  The 
test  requires  that  any  of  the  three  component  equations  should  yield  an  identical  value  for  C. 
This  test  was  also  applied  to  all  the  test  data;  the  results  are  shown  in  Tablel5.  Again,  the 
accuracy  achieved  is  satisfactory.  Table  D-2  contains  details  of  this  calculation. 

The  two  validity  tests  were  successful,  thus  permitting  the  writing  of  the  prediction  equations. 
Example  calculations  for  formulating  prediction  equations  are  shown  in  Section  XIV, 

Volume  11  of  ASD-TR-77-6.  A summary  of  all  prediction  equations  is  listed  in  Table  16. 
Equation  40  is  a combination  of  Eq.  3.4,5  and  corresponding  C.  Equation  41  is  a combina- 
tion of  Eq.  3,6,7  and  corresponding  C,  and  so  on. 


Tabic  15.— Test  of  Validity  for  Constant  Term 


Airplane 

Model 

C = 

SB 

■ 

Average 

Deviation 

(*) 

Value 

Component  Equation  Used 

of 

*2 

TT  wc 

1 VS  2 

1 VS  *3 

TT  vs  TT 

1 VS  4 

B-52 

.1675 

.1631 

.1660 

.1630 

1.59 

.1622 

.1609 

.1651 

.1636 

-.53 

.03394 

.03395 

.0356 

.0338 

2.07 

0.2 

.01658 

.01592 

.01689 

.01591 

3.70 

0.1 

KC-135 

msam 

.7641 

.7470 

.7654 

mm 

mmm 

.3489 

.3485 

.3539 

worn 

0.4 

.0715 

.0694 

.0707 

0.99 

0.2 

.01 504 

.01 549 

.01553 

.01556 

-1.33 

0.1 

F-lll 

.5176 

.5023 

mm 

.5059 

0.34 

0.6 

.2051 

.2077 

H 

.2070 

0.00 

0.4 

.03812 

.03583 

.03553 

3.29 

0.2 

.02085 

.02216 

H 

.02231 

-2.46 

0.15 

70 


l 
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Value  used  to  derive  the  equation 


SECTION  X 

MODEL-TO-SIMULATOR  CORRELATION 


The  prediction  equations  were  next  used  to  correlate  back  witli  the  stopping  distance  data 
collected  in  the  Task  HA  simulation.  A summary  of  errors  in  correlation  is  listed  in  Table 
17.  The  model  to  simulator  correlation  tables  foi  all  /a  conditions  are  included  in  Section 
XIV,  Volume  11,  ASD-TR-77-6. 

The  limitations  (range  of  validity)  of  the  prediction  equations  are: 

• Equations  40,  41 , 45,  46,  49,  and  50  are  applicable  for  p values  of  0.3  to  0.6. 

• Equation  42  is  applicable  for  p range  of  0.2  to  0.3  only. 

• Equations  43,  and  44  are  applicable  for  p range  of  0.05  to  0.2  only. 

• Equations  47,  and  48  are  applicable  for  p range  of  0.1  to  0.3  only. 

• Equations  51 , and  52  are  applicable  for  p range  of  0.15  to  0.3  oniy. 

Fora  given  airplane  model,  the  prediction  equations  are  interchangeable  (alternate  solutions) 
if  their  range  of  applicability  and  validity  is  common.  Thus,  Eq.  40  and  41  are  interchange- 
able, so  are  43  with  44,  45  with  46  and  49  v/ith  50.  Eq.  42  is  a unique  solution  and  not 
interchangeable  with  any  of  its  counterparts,  Eqs.  40,  41 , 43  or  44.  Equat  .ons  47  and  48 
should  have  been  interchangeable;  however,  the  wide-variation  in  the  simulator  baseline  data 
(as  explained  in  Section  IV)  prevented  the  achievement  of  a ± 5 % correlation  accuracy,  i.e., 
ability  of  one  equation  to  correlate  with  the  data  of  the  counterpart  equation.  By  same 
reasoning  Eq.  51  and  52  should  have  been,  but  are  not,  interchangeable. 

Some  airplane  systems  need  only  one  prediction  equation  to  define  the  entire  range  of  p 
values  tested  on  'lie  simulator;  others  needed  more  than  one  equation.  The  reason  for  this 
can  be  comprehended  by  studying  braking  distance  efficiency  curves  for  the  various  systems 
as  shown  by  Figure  24. 

Braking  distance  efficiency,  ns,  is  defined  as  the  ra?:o.  ■ f the  perfect  braking  distance  to  the 
braked  airplane  dictan  e resulting  from  the  simulation. 

ns  = Xp/Xax  100% 

where  ns=  braking  distance  efficiency 

Xp  = perfect  braking  distance 
Xa  = airplane  braking  distance 

The  perlect  braking  distance  is  the  distance  required  to  stop  the  airplane  if  it  is  brakeu  for 
the  entire  jtop  with  maximum  available  braking  force..  Braking  distance  efficiency  indicates 
the  degree  to  which  the  system  meets  its  priina  y requirement  of  stopping  the  aircraft. 


Table  17—  Summary  of  Percentage  Errors 


-4.43  to  +4.91 


-0.77  to  +2.08 
-4.71  to  +1 .91 

-3.00  to  +1 .66 
-1 .22  to  +3.49 

-4.85  to  +4.86 
-9.95  to  +0.00 


-4.58  to  +3.58 
-0.50  to  +9.56 


The  skid  control  systems  for  the  three  subject  airplanes  encompass  two  enerations  of  tech- 
nology namely,  the  old,  and  an  intermediate  type  (see  reference  1 page  93),  The  B-52, 
Hydro-Aire  MK  1 system  was  developed  under  the  old  technology,  the  KC-1 35,  Hydro-Aire 
MK  11  System  from  the  intermediate  technology  while  the  F-l  1 1 . Goodyear  System  is  very 
similar  to  the  Hydro-Aire  MK  11  system  in  design.  Two  well  defined  discontinuities  occur 
in  the  curve  for  the  B-52  for  reasons  related  to  torque-limited-braking,  weight  and  velocity 
combination  and  the  resulting  skid  activity  from  being  in  a particular  region  of  mu-slip  curve 
(see  Section  IV  of  this  report)  while  sharp  changes  in  slope  appear  in  the  curves  for  the  other 
two  systems  at  p values  of  0.3.  Piecemeal  linearization  is  required  when  writing  mathemati- 
cal relationships  for  curves  of  this  nature.  This  in  turn  leads  to  several/multiple  prediction 
equations  for  each  of  the  curves  shown. 

The  sudden  drop  in  the  efficiency  at  0.3p  for  the  B-52  airplane  is  the  result  of  torque  limited 
braking  experienced  at  mu  levels  of  0.3  and  higher;  the  torque  limiting  being  the  result  of  an 
under  designed  brake.  Under  torque  limiting  conditions  there  is  no  antiskid  cycling  because 
the  brake  does  not  have  the  capability  to  counteract  the  maximum  available  ground  force 
fully.  The  prediction  equation  then  becomes  independent  of  the  nj  (i.e.,  mu)  term  and 
equation  2a  takes  the  form: 


(irj)  = F (»3,  rr4) 

(2b) 

or  (ffi)  = K(ir3)a(jr4)P 

The  data  generated  on  the  simulator  for  cases  where  torque  limiting  was  experienced  could 
be  handled  by  this  simplified  form  of  prediction  equation;  however,  the  general  prediction 
equations  provided  in  the  report  for  B-52  would  also  handle  these  cases.  Torque  limiting 
was  experienced  only  at  very  high  weight  and  high  friction  (/i)  value  combinations  for  both 
the  KC-135  and  the  F-l  1 1 airplanes.  The  torque  limiting  cases  for  the  KC-135,  and  F-l  1 1 
were  therefore  excluded  from  the  formulation  of  equations  and  the  correlation. 

The  correlation  data  error  summary  (Tablel7)indicates  that,  for  almost  al!  conditions,  a 
predict; in  accuracy  of  ± 5%  can  be  achieved. 

Even  though  in  the  correlation  process,  comparison  was  made  between  predicted  and  actual 
7r  j values,  that  is: 

(sg/y2)pre(j  versus  (sg/y 2)acjUai 

it  is  tantamount  to  comparing  the  braking  stop  distances  since  both  terms  use  identical  g and 
v values  and  the  distance  term  has  no  exponent. 


SECTION  XI 

ADDITIONAL  CORRELATION  AND  WET  RUNWAY  ANALYSIS 


\ 


l 


5 


1 . PREDICTED  MU  FROM  FLIGHT  TEST  DATA 

The  credibility  of  the  simulator  procedure  was  established  by  comparing  the  simulator  and 
airplane  flight  test  data  and  showing  that  similar  trends  were  obtained  under  identical  condi- 
tions (see  Section  VI).  The  credibility  of  the  prediction  model  has  been  established  by  ob- 
taining a ± 5 correlation  accuracy  m predicting  simulator  stopping  distances.  The  next 
logical  step  is  to  determine  if  the  airplane  flight  test  data  could  be  correlated  to  the  predic- 
tion model.  This  three  way  correlation  process  is  depicted  graphically  m figure  25.  T,ie 
results  of  this  exercise  are  shown  in  Tablel8,  From  the  type  of  information  available  on  the 
flight  test  data  (Table  1 8)the  only  parameter  that  could  be  calculated  by  the  prediction 
equation  was  the  friction  coefficient.  The  predicted  values,  for  both  dry  and  wet  runway 
conditions,  compare  rather  well. 

2.  WET-RUNWAY  ANALYSIS 

During  Task  I1A  simulation  testing,  a wet  runway  was  simulated  so  that  the  available  ground 
mu  was  programmed  to  vary  with  speed,  (See  Figure  26).  The  mu  values  (end  points)  used 
were  0.05  at  brake  application  speed  and  0.5  at  the  end  of  the  stop.  Additional  wet  runways 
were  also  sii  dated  with  end  points  being  0.05  to  0.4  and  0.05  to  0.3. 

The  average  .-.tie  of  peak  available  mu  for  the  braking  system  was  tin’  lown,  so  it  was 
decided  to  use  the  component  equations  formed  eariier  (the  w j veisus  7r->  relationships)  to 
calculate  peak  available  mu.  Based  on  calculations  for  wet  runways,  prediction  equations  were 
generated  for  each  wet  runway.  With  these  prediction  equations,  a correlation  prediction 
accuracy  analysis  was  conducted  as  before  and  satisfactory  results  were  obtained.  The  com- 
ponent equations,  the  prediction  equations,  and  the  correlation  error  tables  for  the  wet  run- 
ways tested,  are  included  in  this  section,  (see  Tables  1 9 through  23).  The  details  of  this 
analysis  are  repo  ted  in  ASD-TR-77-6.  Volume  II.  Section  XV.  The  results  show  that  the 
wet  runway  data  based  prediction  equations  relate  with  each  other  with  the  same  accuracy 
as  the  fixed  mu  prediction  equations,  and  give  additional  confidence  to  the  selected  method- 
ology for  forming  prediction  equations. 
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Config.  & peri. 


Flight 

Test  compare: 

Braking  distance 
Skidding  pressure 
Number  of  skids 
Depth  of  skids 
Rate  in  and  out  of 
skids 

General  control 


Comparison 


Simulator, 
Duplicate  a known 
flight  test  cond. 

• A/C  weight 

• C.G.  location 

• Brakes  on  speed 


Table  18—Comparisc.,)  with  Flight  Tesi  Data 


VALUE  OF  PEAK  AVAILABLE  MU  USED  ON  SIMULATOR  FOR  CORRELATION  WITH  FLIGHT  TEST  DATA 
CALCULATED  MU  =AvG./nBRA|aHG 

WHERE^avg.  IS  THE  VALUE  TAKEN  FROM  TABLE  I,  REFERENCE  { 15). 

AND  nD0.„TMr  IS  THE  BRAKING  DISTANCE  EFFICIENCY  DETERMINED  BY  SENSITIVITY  STUDY 
DKAKInG 


PERCENT  BRAKE  APPLICATION  SPEED  V,lo  BRAKE 


Table  19— Summary  of  Component  Equations  ■ Wet  Runways 


Equation 


, , x -2.2504 


= .1  747  (tt2) 

= 4.8081  (tt  3) 

= 24.5495  (tt4) 


.1  3006- . 07741 %SP 
v 17.75AO  - .1836 


= 2.3524  (it 2) 

= 5.8154  (tt  ) 

= 32.7278  (ir4) 


v . 1 275- . 1 258%SP 


19.6CAp  -.1991 


.05  to  .3  K)  = 2.3524  (n  ) 
(.127)  1 Z 


-.5188 


(tt  ) = 6.8597  (tt3)  . 1 049- . 1 24 3%SP 


(*,)  = 37.7863  (tt4) 


20.97&P  - .19788 


KC-135  .05  to  0.5  (tt,  ) = .5893  K) 

(.166;  ' 6 

(tt,  ) = 4.4163  (tf3) 

(^}  = 15.2166  (ir4J 


-1.1410 


.1994-.1 225%SP 
lS.OOAp  - .1063 


.05  to  0.4  (it,)  = .5893  (^J  '1J410 

(''54)  (,,1  - 4.8494  (,3) 

(tt,  ) = 11  .1128  (tt4)  21  .OOAp  - .07158 

.05  to  0.3  (tt,  ) = .5893  (t^)  "1,141° 

('137)  1 \ R c7Cji  in  ^ -1913-.1731  %SP 

(tt,  ) = 5.67S4  ( tt3) 

t \ ,c  ,40,  , \ 17. 54Ap -.08797 

K)  = 15.3497  (tt.) 


* was  calculated  with  £q.  (11)  data  plus  wet  runway  data. 


Table  19— Summary  of  Component  Equations  - Wet  Runways  (Concluded) 


Ai rplane 
model 


Equation 


F-lll 

.OS  to  .5 

(*,)  = .9467  (n2)  ‘1’0490 

(IT,)  = 4.3767  (1T3)  -1  790-.0834%SP 

(ir,)  = 33.7671  (tt4)  8*0aP-'1606 

(31) 

(.231) 

(63) 

(64) 

.05  to  .4 
(.218) 

(it,)  = .9467  (tt2)  -1-0490 

«,)  - 4.6220  (.3)  •1976-1°46^P 

(it,)  = 46.5388  (*4)  7-75AP--1804 

(31) 

(65) 

(66) 

.05  to  .3 

(ir,)  = .9467  (n2)  -1-0490 

(it,)  = 5.2915  (tt3)  .1973  - . 1 296%SP 

(tt,)  = 55.0660  (ir4)  ‘-1827 

(31) 

(.193) 

(67) 

(68) 

Table  21—  Test  of  Validity  for  Constant  Term  - Wet  Runways 


Airplane 

Model 

C = 

O'f'Y  ’3- 

v]2 

■ 

Average 

Deviation 

(%) 

Value 

of 

"2 

Component  Equation  Used 

*1  vs  *2 

*1  VS  ^3 

1T1  vs  *4 

B-52 

.04395 

.04349 

.04362 

.04373 

0.10 

.023 

.02954 

.02957 

.03112 

.02954 

1.83 

.175 

.02124 

.02121 

.02373 

.02129 

3.62 

.127 

KC-1 35 

' 

■i 

.04906 

— 

.166 

■SB 

.04146 

■BSH 

1 .10 

.154 

..03065 
— 1 

m 

.031 56 

■ 

0.99 

.137 

F-ll  1 

mm 

.05224 

mm 

0.99 
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Table  23.— Summary  of  Percentage  Errors  - Wet  Runway: 


ff 


AIRPLANE 

USING  PRED. 

APPLIED  TO 

ERROR 

MODEL 

EQ.  FOR  n2 

DATA  0 ^ 

RANGE 

(EQ.  NO.) 

B-52 

0.23 

0.23 

-2.98  to  +3.58 

(42) 

0.175 

0.175 

-4.25  to  +4.05 

(69) 

0.127 

-7.08  to  +2.64 

0.127 

0.127 

-7.95  to  +2.76 

(70) 

0.175 

-5.40  to  +4.69 

KC-1 35 

0.166 

!■ 

-4.76  to  +4.73 

(71) 

mXEM 

-9.88  to  +7.30 

0.137 

-5.30  to  +6.00 

0.154 

0.154 

-4.04  to  +5.58 

(72) 

0.166 

-5.68  to  +6.47 

0.137 

-6.46  to  +5.69 

0.137 

0.137 

-5.32  to  +5.13 

(73) 

0.166 

-7.46  yo  +5.91 

0.154 

-6.24  to  +5.37 

F-lll 

0.231 

0.231 

-4.68  to  +3.67 

(74) 

0.218 

-1.93  to  +4.34 

0.193 

-0.75  to  +b.84 

0.218 

0.213 

-2.21  to  +4.69 

(75) 

0.231 

-5.50  to  +3  79 

0.193 

-4.43  to  +3.17 

0.193 

0.193 

-3.92  to  +2  98 

(76) 

0.231 

-6.13  to  +3.38 

0.218 

-3.15  to  +1 .34 

Table  23.— Summary  of  Percentage  Errors  - Wet  Runways 
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SECTION  XU 


TBPS  COMPLETENESS  AND  VERIFICATION 

1.  COMPATIBILITY  VERIFICATION 

Under  Task  I analysis  the  Friction  Prediction  Subsystem  concept  was  developed  and  is 
shown  here  in  figure  27,  Validation  of  the  Braking  Performance  Prediction  Model  was 
established  under  Task  11  and  the  resulting  Braking  Prediction  Subsystem  concept  is 
depicted  in  figure  28.  These  two  subsystems,  when  integrated,  will  form  the  Total  Braking 
Prediction  System  as  shown  in  figure  29.  The  compatibility  between  the  output  of  the 
FPSS  and  the  input  to  the  BPSS  is  solely  dependent  upon  the  tire-model  verification  which 
in  turn  is  based  upon  conducting  the  recommended  tire  testing.  However,  based  on  our 
engineering  judgment  and  available  tire  test  data,  there  is  no  reason  to  believe  that  this 
compatibility  will  not  be  achieved. 

Assuming  that  the  recommended  tire  test  data  will  be  obtained  and  that  the  output  of  the 
FPSS  will  be  in  the  form  of  a family  of  mu-velocity  curves,  we  can  proceed  to  the  TBPS 
analysis  of  the  simulator  conditions.  As  explained  in  section  XI,  wet  runway  analysis,  the 
wet  runway  curve  was  input  to  the  simulator  as  a mu-velocity  curve  (figure  26).  In  accord- 
ance with  the  suggested  TBPS  and  its  subsystems,  three  wet  runway  curves  (.05  to  .5,  .05 
to  .4,  .05  to  .3)  forming  the  so-called  family  of  curves  were  input  for  each  of  the  three  air- 
planes studied.  The  sensitivity  tests  were  run  and  the  data  analyzed,  as  already  explained  in 
section  XI.  The  prediction  equations  obtained  for  each  of  these  wet  runways  and  their 
cross-correlation  accuracy  were  shown  to  be  in  conformity  with  the  fixed  mu  equations  and 
correlation  for  each  of  the  respective  airplanes.  This  clearly  establishes  that  when  the 
required  significant  parameters  are  input  into  the  BPSS,  the  correct  stopping  prediction 
will  result. 


2.  FINALIZATION  OF  TBPS  SPECIFICATION 

The  specification  for  the  TBPS  is  comprised  of  the  FPSS  performance  specification,  the 
general  BPSS  model  and  the  tire  correlation  model.  The  basic  concept  of  the  FPSS  is  shown 
in  Figure  27,  The  FPSS  performance  specification  was  established  under  Task  I of  this 
contract  and  discussed  in  detail  in  ASD-TR-77-7,  The  FPSS  specification  also  covered  the 
system  fabrication  criteria.  The  general  BPSS  model  was  established  and  validated  by  a 
sensitivity  study  of  eight  different  airplane  braking  systems.  The  general  concept  of  the 
BPSS  is  depicted  in  Figure  28.  Certain  hardware  and  data  are  required  for  a successful 
simulation  of  any  given  aircraft  braking  system. 

Basic  vehicle  system  parameters  and  brake  system  hardware  are  incorporated  into  the 
computer/hardware  simulation  for  each  of  the  airplanes  being  evaluated.  The  basic  data 
are  available  from  existing  specifications,  qualification  reports,  and  flight  handbooks,  and 
are  generally  no  cost  items. 
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Figure  27. —Friction  Prediction  Subsystem  ( FPSS ) 


Figure  28—  Braking  Prediction  Subsystem  (BPSS) 


Brake  control 
simulator 


The  requirements  can  be  divided  into  three  categories  as  follows: 

• Basic  vehicle  parameters 

• Landing  gear  strut  parameters 

• Tire,  wheel  and  brake  system  parameters 

The  specific  requirements  for  each  category  are  listed  in  Table  24.  Any  of  the  less  critical 
parameters  which  are  not  available  can  be  estimated  for  the  simulation. 

Brake  system  hardware  required  for  the  USAF-F-4,  C-141 , B-52,  KC-135,  and  F-l  1 1 simula- 
tions are  listed  in  Table  25.  Items  unique  to  each  skid  control  system  are  checked  accordingly. 

In  addition  to  the  above  requirements,  airplane  flight  test  records  showing  brake  pressure, 
wheel  speed  valve  voltages  and  strut  loads  are  very  helpful  and  necessary  for  establishing 
the  credibility  of  the  individual  simulation. 

The  tire  correlation  model  was  also  established  under  Task  I analysis  and  ihe  concept  is 
shown  here  in  figure  30.  The  integration  of  the  three  subsystem  concepts  has  already  been 
discussed  and  shown  in  figure  29.  The  information  provided  herein  and  in  ASD-TR-77-7, 
together  should  be  sufficient  to  allow  total  system  fabrication  and  development. 

3.  VERIFICATION  TEST  PROGRAM 

An  extensive  test  program  must  be  conducted  for  the  verification  of  the  Total  Braking 
Prediction  System.  It  involves  the  use  of  the  Friction  Prediction  Subsystem  (ground  vehicle), 
the  Prediction  Equation,  Tire  Correlation  Method,  and  at  least  one  suitable  aircraft.  The 
overall  goal  is,  of  course  the  prediction  of  braking  distance  with  the  methods  and  tools 
recommended  during  the  course  of  this  program.  Comparison  of  predicted  braking  distance 
to  the  value  measured  from  actual  airplane  tests  will  give  a measure  of  correlation.  Measure- 
ments taken  with  the  ground  vehicle  and  actual  airplane  braking  test  must  be  conducted 
within  a few  minutes  of  each  other  so  that  runway  and  environmental  conditions  can  vary 
the  least  possible  amount.  The  major  tasks  which  must  be  accomplished  are  shown  on  Figure 
31 . A ground  vehicle  which  is  configured  and  performs  according  to  the  Friction  Prediction 
Subsystem  specification  is  used  to  conduct  friction  evaluation  tests  on  the  runway.  With  the 
use  of  the  previously  developed  prediction  equation  for  the  vehicle  a friction  coefficient  can 
be  determined.  Due  to  the  physical  and  operational  differences  between  vehicle  and  airplane 
tires  a correction  factor  must  be  applied  to  the  friction  coefficient  calculated  from  the 
ground  vehicle  test  before  it  can  be  applied  to  the  airplane  braking  distance  prediction 
equation.  The  airplane  conditions  at  brake  application  are  to  be  used  in  the  prediction  equa- 
tion. Hence  in  this  verification  program  a braking  distance  comparison  can  only  be  made 
after  the  airplane  and  ground  vehicle  tests  are  performed.  In  actual  operational  use  the 
anticipated  condition  at  touchdown  will  be  used  for  a prediction  of  braking  distance.  The 
sensitive  parameters  such  as  brake  application  speed  and  friction  coefficients  can  then  be 
varied  to  determine  the  margin  of  the  runway. 
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Table  24.— Basic  Airplane  Data  Required 


• Airplane  parameters 

1.  Location  of  the  landing  gear  relative  to  the  center  of  gravity. 

2.  Airplane  landing  gross  weight  range,  (including  corresponding  stall  speeds). 

3.  Airplane  center  of  gravity  range  (vertical  and  horizontal)., 

4.  Effective  wing  area. 

5.  Aerodynamic  coefficients— lift  and  drag  (nominal  landing  flaps,  spoilers,  up  and  down). 

6.  Total  engine  idle  thrust  versus  speed. 

7.  Airplane  mass  moment  of  inertia  about  the  center  of  gravity  in  the  pitch  direction. 

• T ire,  wheel  and  brake  system  parameters 

1.  Weight 

a)  Main  gear  tire  and  wheel  assembly. 

b)  Total  brake  assembly. 

c)  Brake  heat  sink. 

2.  Brake  mean  torque  versus  pressure. 

3.  Number  of  braked  wheels. 

4.  T<re  size  and  ply  rating. 

5.  Tire  mechanical  properties. 

6.  Brake  hydraulic  system  diagrams  showing  line  sizes,  line  lengths,  and  materials. 

7.  Mass  moment  of  inertia  of  the  wheel,  brake  and  tire  assembly  about  the  axle  centerline. 

• Landing  gear  strut  parameters 

1 ..  Vertical  spring  rate  of  the  main  gear  oleo. 

2.  Vertical  spring  rate  of  the  nose  gear  oleo. 

3.  Effective  length  of  the  main  gear  strut. 

4.  Fore-aft  spring  rate  of  the  strut. 

5.  Strut  fore-aft  natural  frequency  range. 

6.  Total  effective  mass  of  the  strut. 

7.  Main  landing  gear  layout  drawing  showing  basic  dimensions. 

8.  Truck  size,  weight,  center  of  gravity  and  pitch  moment  of  inertia  (where  applicable). 

9.  Fore-aft  mass  moment  of  inertia  of  main  gear  strut. 
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BRAKE  SWIVEL  ASSEMBLY 


Figure  30— Tire  Correlation  Model 


Vehicle  stopping  tests  on  runway  from 
70  mph  to  full  stop  to  evaluate  friction 
coefficient  of  runway.  Performed  prior 
and  after  airplane  test 


A non-dimensional  equation  based  on 
variables  of  the  ground  vehicle  developed 
specifically  to  calculate  peak  friction 
coefficient 


Tire  size,  speed-range,  and  inflation 
pressure  of  ground  vehicle  ana  airplane 
are  not  similar.  A correlation  factor, 
which  must  yet  be  developed  according 
to  a correlation  concept  discussed 
earlier 


Applying  the  correlation  factor  to  modify 
the  friction  coefficient  derived  from  the 
ground  vehicle  test 


Airplane  braking  prediction  equation 
utilizing  the  modified  friction 
coefficient  and  airplane  parameters 
which  are  recorded  at  touchdown  during 
the  airplane  test 


For  this  method  to  be  successful  the 
calculated  braking  distance  should  be 
within  5%  of  the  measured  distance 


The  distance  recorded  with  nose  wheel 
revolution  counter  from  Drake  applica- 
tion to  full  stop 


Figure  31.— Block  Diagram  for  Verification  Test  Program 


TEST  HARDWARE 


The  following  major  items  are  required  for  conducting  the  test  program: 

a.  Friction  Prediction  Subsystem -(ground  vehicle) 

A ground  vehicle  with  performance  as  described  in  the  specification  (Task  I Report,  ASD- 
TR-77-7)  must  be  available  for  the  measurement  of  ground  friction.;  The  most  relevant 
features  of  the  ground  vehicle  are  that  the  test  tire  is  an  aircraft  tire  which  can  be  braked 
with  an  antiskid  system.  Performance  of  the  antiskid  system  must  be  similar  to  that  of  an 
advanced  type  aircraft  antiskid  system.  1 he  vehicle  is  to  be  used  for  an  assessment  of  runway 
friction  which  will  be  used  for  the  prediction  of  aircraft  stopping  distance.  A prediction- 
equation-Mu  is  part  of  the  Friction  Prediction  Subsystem  which  allows  the  calculation  of 
Mu  based  on  the  brake  application  speed  and  the  resulting  stopping  distance..  Other  variables 
such  as  environmental  conditions  may  also  be  involved. 

b.  Airplane  and  Prediction  Equation 

A suitable  aircraft  equipped  with  an  antiskid  system  must  be  made  availabl"  and  perform 
braking  stops  for  which  the  distance  will  be  predicted  from  FPSS  measurements  and  the 
Prediction  Equation.  It  is  recommended  that  the  Advanced  737  (T-43A)  be  used  for  this 
program  for  the  following  reasons: 

• Proven  short  field  capability 

• Excellent  braking  performance 

• Efficient  modern  antiskid  system 

• Economic  test  airplane  due  to  two  engine  configuration 

• Braking  Prediction  Equation  applicable  over  a wide  range  of  friction  coefficier.'0 

The  T-43A  has  proven  short  field  capability  and  tests  could  be  conducted  at  a variety  of 
different  airfields.  Should  an  extremely  long  runway  be  available,  the  braking  tests  can  be 
run  in  an  accelerate  stop  fashion  without  the  aircraft  ever  becoming  airborne.  Special  brake 
cooling  sources  would  have  to  be  provided  because  repeated  braking  would  result  in 
overheating  of  the  brakes. 

The  T-43A  is  equipped  with  a modern,  fully  modulating  antiskid  system  with  efficient  opera- 
tion over  all  ranges  of  operating  conditions.  The  operation  of  the  antiskid  system  is  repre- 
sentative of  that  realized  on  the  latest  and  most  modern  jet  aircraft. 

The  T-43A  is  an  economic  test  bed  due  to  the  two  engine  design.  Fuel  and  maintenance 
expense  will  be  less  than  those  encountered  on  three  and  four  engine  aircraft. 

During  the  course  of  the  Combat  Traction  II,  Phase  II  Program,  a Braking  Prediction  Equa- 
tion has  been  developed  for  the  T-43A  which  covers  a wide  range  of  friction  coefficients. 
Good  correlation  to  simulator  results  and  flight  test  data  has  been  demonstrated. 
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c.  Tire-Correlation  Method 


The  test  tire  on  the  ground  vehicle  is  similar  in  design  and  applied  antiskid  operation  to  the 
aircraft  tire.  However,  the  size,  inflation  pressure,  and  speed  range  over  which  the  tires  are 
being  used  are  different.  It  is  anticipated,  therefore,  that  the  friction  coefficient  derived 
from  the  ground  vehicle  must  be  modified  before  being  used  in  the  airplane  braking  distance 
prediction  equation.  Based  on  a number  of  isolated  data  points  for  different  aircraft  tires, 
a correlation  concept  has  been  developed.  However,  a much  wider  data  base  must  be 
developed  before  this  correlation  concept  can  be  used  with  any  degree  of  confidence.  This 
tire  correlatibn  method  must  be  developed  before  the  Total  Braking  Prediction  System  veri 
fication  test  program  is  initiated. 

-»>> 

TEST  PLANE 

The  major  portion  of  the  test  consists  in  making  performance  stops  on  both  wet  and  dry 
runways  for  a variety  of  initial  airplane  conditions.  These  include,  but  are  not  limited  to, 
the  conditions  shown  on  the  matrix  of  figure  32.  The  spoiler  UP  and  DOWN  configuration 
will  result  in  the  widest  possible  variation  on  aerodynamic  lift  and  drag.  The  overspeed 
touchdown  conditions  are  very  typical  of  the  operational  environment  airplanes  encounter 
every  day.: 

Ideally,  testing  should  be  conducted  under  calm  wind  conditions,  however,  as  an  upper 
limit  a steady  wind  of  no  more  than  5 knots  from  any  direction  is  acceptable.  In  flying  the 
airplane  to  touchdown  a normal  (2.5  degree)  glide  slope  should  be  maintained,  touchdown 
should  be  aimed  at  the  1000  ft  marker.  During  the  flare  the  pilot  should  cut  the  engine  power 
to  idle.  At  touchdown  the  spoilers  (speed  brakes)  should  be  extended  (when  required)  as  the 
nose  gear  is  lowered  to  the  ground.  Then  the  brakes  should  be  applied  smoothly  and  firmly 
until  maximum  metered  brake  pressure  is  reached  and  maintained  to  a full  stop.  Tc  limit 
the  transition  distance  the  time  elapsed  from  touchdown  to  brake  application  should  not 
exceed  four  seconds.  No  reverse  thrust  should  be  applied.  This  technique  shall  be  used  for 
both  wet  and  dry  runway  tests. 

The  aircraft  shall  be  instrumented  for  the  recording  of  at  least  the  following  parameters: 


a.: 

Wheel  speed  transducer  signal 

(four  locations) 

b. 

Brake  pressure  and  metered  pressure 

(eiftht  locations; 

c. 

Antiskid  valve  voltage 

(fout  locations) 

d. 

Inboard  brake  center  stator  temperatures 

(twc  -oca.  ions) 

e. 

Nose  gear  revolution  counter 

(o>;e  location) 

f. 

Engine  pressure  ratio 

(t  vo  locations) 

g- 

Spoiler  handle  position 

(two  locations) 

h. 

Braking  distance  -Theodolite 
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Floors  32. -Test  Matrix 


For  the  dry  runway  test  no  special  runway  preparation  is  required  as  long  as  there  are  no 
excessive  rubber  or  dirt  deposits.  Considerably  more  effort  is  required  to  prepare  the  “wet” 
runway  where  no  rainfall  is  available.  When  artificial  wetting  of  the  runway  is  required  this 
is  best  accomplished  with  watertrucks,  with  a capacity  of  5,600  gallons  per  truck.  One  hour 
prior  to  each  series  of  wet  runway  tests  the  runway  must  be  prewetted  with  14  water  trucks. 

As  the  truck  proceeds  down  the  runway,  the  water  can  be  dumped  on  the  runway  surface 
through  one  or  two  nozzles,  on  each  truck,  flooding  the  runway  briefly  before  it  runs  - ill.. 

The  actual  water  discharge  rate  and  truck  speed  depends  to  some  extent  on  the  runway 
surface  textpre  and  must  be  adjusted  as  required.  Before  each  test  additional  6 or  7 trucks 
dump  water  to  either  side  of  the  runway  centerline.  The  test  conductor  must  ensure  that 
the  proper  length  of  the  runway  is  covered  with  water.  Because  the  water  will  gradually 
run  off  the  runway  readings  with  the  ground  vehicle  shall  be  taken  just  before  and  immediately 
after  the  airplane  stopping  test.  The  ground  vehicle  test  consists  of  accelerating  to  the 
desired  speed  and  then  braking  to  a full  stop  with  the  antiskid  braked  tire.  This  test  must  be 
repeated  for  the  length  of  the  test  runway.  Runway  friction  values  can  then  be  derived 
from  the  prediction  equation  and  data  recorded  during  the  vehicle  test.  The  friction  values 
shall  be  averaged  when  they  are  applied  to  the  prediction  equation  of  the  airplane  to  com- 
pensate for  the  actual  time  at  which  the  airplane  touched  down.  Close  timing  of  these 
readings  is  essential  to  make  the  results  representative  of  those  that  the  airplane  encounters. 
Communication  by  radio  between  the  test  coordinator,  pilot,  ground  vehicle  driver  and  the 
water  trucks  convoy  is  essential.  Air  temperature,  wind  velocity  and  direction  shall  be  moni- 
tored at  two  points  on  the  side  of  the  runway,  approximately  1 2 ft  above  the  runway  surface. 
One  monitoring  point  should  be  located  near  the  touchdown  point,  the  other  about  2000  ft 
down  the  runway. 

TEST  CONDITIONS 

The  recommended  test  conditions  are  shown  on  Figure  32  and  include  variations  in  airplane 
touchdown  speeds,  weight,  aerodynamic  configuration  with  spoiler  UP  and  DOWN,  and  dry 
and  wet  runway.  A total  of  31  tests  are  recommended.  Due  to  the  resulting  long  braking 
distances  caused  by  touchdown  with  excessive  speed  and  undeployed  spoilers  on  wet  run- 
ways, some  of  those  conditions  are  not  recommended.  If  touchdown  is  accomplished  at  the 
1000  ft  marker  and  1000  ft  is  consumed  in  derotating  the  airplane  and  applying  brakes,  then 
more  than  1000  ft  of  margin  will  remain  for  the  longest  stop  as  calculated  with  the  737 
prediction  equation  when  testing  is  conducted  on  a 10,000  ft  wet  runway. 

CORRELATION  CRITERIA 

For  this  evaluation  test  the  braking  distance  can  only  be  calculated  after  the  airplane  test 
because  the  airplane  brake  application  speed  can  only  be  determined  after  the  test  records 
are  examined.  Also  for  wet  runway  tests  the  airplane  friction  coefficient  will  be  available 
after  the  second  test  which  is  conducted  shortly  after  the  airplane  test.  The  method  of 
correlation  is  deemed  successful  if  calculated  distances  are  within  5%  of  those  measured 
with  the  airplane  tests. 

In  operational  use  the  anticipated  brake  application  speed  will  be  used  based  on  the  touch- 
down weight  of  the  airplane  and  the  friction  coefficient  will  be  derived  from  a calibration 
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chart  based  on  amount  of  water  on  the  runway.  The  friction  calibration  chart  must  be  veri- 
fied or  corrected  periodically  by  comparing  it  to  the  results  obtained  from  ground  vehicle 
tests. 

TEST  EXTENTION 

The  preceding  test  plan  was  formulated  specifically  for  the  T-43A  (Advanced  737).  However, 
it  can  be  expanded  to  include  other  aircraft.  Additional  testing  with  satisfactory  agreement 
between  measured  and  calculated  braking  distance  can  only  help  to  increase  the  level  of 
confidence  on  theTotal  Braking  Prediction  System.  Additional  aircraft  are  available  for  which 
a Braking  Prediction  Equation  has  been  developed.  These  are  the  727,  747,  KC-135,  B-52, 
C-141 . F-l  1 1 and  F-4.  It  is  anticipated  that  the  best  agreement  between  predicted  and 
measured  braking  distance  will  exist  for  those  airplanes  equipped  with  an  advanced  antiskid 
system  which  operates  efficiently  over  a wide  range  of  conditions.  In  planning  tests  with  the 
airplanes  it  must  be  recognized  that  the  range  of  variables  such  as  airplane  touchdown  speed 
and  runway  conditions  are  limited  by  the  runway  length.  The  airplanes  listed  above  land  at 
higher  speeds  which  result  in  longer  braking  distances  than  the  T-43A.  Hence,  the  range  of 
variables  will  be  more  limited. 


SECTION  XIII 
SUMMARY  OF  RESULTS 


1.  B-52 


Application  of  brakes  at  the  same  speed  regardless  of  the  gross  landing  weight  or  touchdown 
speed  (90  Knots  for  G and  H models,  70  Knots  for  earlier  versions)  is  an  operating  procedure 
that  is  unique  to  the  B-S2.  Absence  of  weight  as  an  independent  variable  in  the  prediction 
model/equation  required  that  a correction  factor  be  calculated  to  convert  a given  weight 
change  into  an  equivalent  velocity  change  (see  Appendix  B)  for  all  mu  levels. 


A correction  factor  was  also  needed  to  convert  a given  wind  component  into  an  equivalent 
velocity  component  to  be  used  in  the  prediction  equation  (Appendix  A)  for  all  mu  levels. 


The  exponent  for  the  *4 
variation  at  all  mu  levels. 


term  had  to  be  modified  to  account  for  air  density 
(Appendix  E). 


The  brake  system  on  the  B-52  operates  in  three  distinct  regions,  a rapid  skid  cycling  region 
(H  = .05  to  .2),  a transition  regie  (ji  = .2  to  .3)  involving  a combination  of  skid  cycling  and 
torque  limiting  and  a complete  torque  limited  operation  (/i>0.3).  The  torque  limiting  is  the 
result  of  an  underdesigned  brake  and  braking  distance  becomes  independent  of  available  mu. 
Hence,  the  prediction  equation  may  be  simplified  to:  (also  see  equation  2b,  Section  X) 


The  numerical  methods  utilized  in  calculating  the  component  equations  and  prediction  equa- 
tions showed  that  for  torque  limited  cases  the  exponent  of  *2  (mu)  approached  a value  of 
zero.  Hence  the  prediction  equation  became  independent  of  mu  and  this  result  coincides 
with  the  physical  nature  of  torque  limiting. 

2.  KC-135 

Lack  of  repeatability  in  baseline  distances  on  the  simulator  was  traced  back  to  variations  in 
anti-skid  valve  characteristics.  Extensive  testing  of  the  antiskid  valve  indicated  that  the 
allowed  tolerance  band  on  the  valve  performance  was  not  compatible  with  the  target 
distance  scatter  of  1 to  2%  on  the  brake  control  simulator.  This  in  turn  prevented  the 
desired  ± 5%  accuracy  in  data  correlation  and  cross-correlation,  especially  at  low  mu’s. 

The  correction  factor  for  a given  wind  component  was  needed  only  at  low  mu’s  (<0.3) 
conditions. 

The  modified  *4  exponent  for  density  variation  was  necessary  only  at  low  mu’s  (<0.3). 
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Torque  limited  braking  was  experienced  in  only  3 out  of  88  test  conditions  (high  wt-high 
velocity  combination)  and  were  excluded  from  the  correlation  process. 

3.  F-lll 

Due  to  lack  of  repeatability  in  baseline  braking  distances  and/or  wheel  lockups  occurring  at 
0. 1 mu  level,  0. 1 5 mu  was  established  as  the  lower  mu  limit  for  the  F-lll.  The  data  scatter 
m the  skid  control  valve  performance  was  not  as  wide  as  for  theKC-135  valve  but  was  enough 
to  substantially  exceed  the  1 to  2 7c  desired  repeatability.  As  a result,  the  desired  ± 5%  cross- 
correlation of  data  was  not  achieved  at  low  mu  levels  (<0.3p). 

The  effect  of  varying  weight  alone  was  inconsistent  as  it  sometimes  increased  and  sometimes 
decreased  the  stopping  distance..  A general  trend  is  however  evident  (Figure  2 1 )., 

The  correction  factor  for  a given  wind  component  was  needed  only  at  low  mu’s  {<0.3). 

The  modified  it 4 exponent  for  density  variation  was  necessary  only  at  low  mu’s  (<0.3). 

Torque  limited  braking  was  experienced  in  only  3 out  of  88  test  conditions  (high-wt,  high- 
velocity  combination)  and  were  excluded  from  the  correlation  process. 

4.  GENERAL 

One  of  the  requirements  of  the  Task  H analysis  was  to  compare  the  results  (prediction  model) 
of  the  current  study  with  those  of  the  earlier  sensitivity  analysis,  reference  ( 1 ),  and  establish 
compatibility  if  necessary  by  modifying  the  earlier  models.  The  idea  was  to  have  one  general 
model  applicable  to  all  study  airplanes.  This  exercise  showed  the  need  for  modifying  the  jr^ 
exponent  for  density  variation  for  F-4  airplane  at  all  mu  levels  and  the  wet  runway  equations 
for  the  other  four  (727,  737,  747,  and  C- 141)  airplanes.  The  general  methodology  has 
however  remained  unchanged  in  th3t  the  number  of  variables,  the  number  of  it  terms,  the 
nature  of  prediction  model  and  the  correlation  accuracy  have  remained  the  same.  Validation 
of  the  previously  developed  prediction  model  has  thus  been  established. 


SECTION  XIV 

CONCLUSIONS  AND  RECOMMENDATIONS 


1.  CONCLUSIONS 

The  following  conclusions  can  be  drawn  from  the  analysis  of  the  available  data  generated 

during  this  program  as  well  as  the  previous  contracted  effort  (reference  1 ). 

V 

• Dimensional  analysis  technique  can  successfully  express  braking  distances  in  the  form 
of  an  equation. 

• The  prediction  model/eqaation  can  predict  braking  distances  within  ± 5%  for  the  air- 
planes studied. 

• Simulation  and  analysis  of  additional  military  aircraft  has  helped  validate  the  previously 
developed  prediction  model  as  well  as  establish  compatibility  of  the  model  for  various 
types  of  airplanes  e.g.,  bombers,  transports,  fighters,  etc. 

• Although  correction  factors  may  be  needed  for  density,  wind  velocity  and  weight  alone 
variation  in  some  airplane  prediction  equations,  the  general  methodology  has  been 
shown  to  be  valid. 

• A better  understanding  of  contributions  made  by  unique  operating  procedures,  gear 
geometry,  type  of  skid  control  system  and  operating  range  has  resulted  from  these 
sensitivity  studies. 

• A Total  Braking  Prediction  System  (TBPS)  concept  has  been  defined  and  a suitable 
airplane  test  program  has  been  outlined  to  verify  the  same. 

2.  RECOMMENDATIONS 

Certain  conclusions  and  recommendations  were  listed  in  Task  I report,  ASD-TR-77-7,  that 

must  be  kept  in  an  overall  perspective.  The  analysis  established  that: 

• Tire  test  data  must  be  collected  under  fully  controlled  conditions  in  order  to  validate/ 
modify  the  tire  correlation  model  established  by  dimensional  analysis, 

• A Friction  Prediction  Subsystem  (FPSS)  specification  criteria  could  be  generated  to 
develop  a suitable  ground  vehicle  that  gives  a meaningful  measurement  of  the  tire- 
runway interface  mu. 

• None  of  the  existing  ground  vehicles  meets  the  FPSS  specification  criteria. 

With  this  background  plus  the  conclusions  of  Tasks  II  and  III  listed  earlier  the  following 

recommendations  should  be  carried  out,  preferably  in  the  order  listed: 
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Tire  test  data  must  be  collected  as  outlined  under  Task  1 and  the  tire  model  validated/ 
modified. 

A suitable  ground  vehicle  must  be  developed  based  on  the  specification  criteria 
established  under  Task  I., 

An  extensive  test  program  must  be  conducted  for  the  verificat  on  of  the  TBPS  as  out- 
lined m this  report. 

In  order  for  the  TBPS  and  its  subsystems  (FPSS  and  BPSS)  to  be  operationally  meaning- 
ful, the  following  areas  of  work  have  to  be  resolved  in  addition  to  the  test  programs 
recommended  above: 

• Classification  of  runways 

• Method  of  measuring/indicating  rainfall  intensity/water  depth 

• Runway  monitoring  system  standardization 

• Enacting  and  enforcing  regulations  regarding  proper  maintenance  of  runways/ 
friction  levels. 


APPENDIX  A 

CALCULATION  OF  WIND  FACTOR 


s 


* 

',v 


% 


9 

5 
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Example  for  B-52: 

H = 0.6;  Baseline  brake  application  velocity,  V = 90  Knots 
Procedure: 

V 

For  a 5 knot  tail  wind  component  (-ve),  at  .6p  the  stopping  distance  is  2021  feet  compared 
to  the  baseline  (no  wind)  stopping  distance  of  1 779  feet  as  shown  in  figure  1 7,  section  VIII. 
At  2021  feet,  the  equivalent  brake  application  velocity  is  seen  to  be  97  knots  or  an  increase 
in  baseline  velocity  of  7 knots.  Similarly  a 10  knot  headwind  component  is  equivalent  of 
76.5  knots  brake  application  velocity  or  a decrease  of  13.5  knots.  This  is  illustrated  in  the 
first  three  columns  (table)  below..  The  remaining  calculations  pertain  to  obtaining  an  average 
factor  F tiiat  reflects  change  in  brake  application  velocity  per  knot  change  in  wind  component 
as  a function  of  friction  value  p. 


Wind 

component 

VW 

(knots) 

Equivalent 

velocity 

change 

Aveq. 

(knots) 

(V  -AvEQ ) 
sVEQ. 

v“  ] 
Lv-vw  J 

(AF/Vyy) 

■ 

E 

l 

-10 

-16 

106 

.0600 

-5 

-7 

97 

.0210 

■SB 

.0051 

.6 

+10 

+13,5 

76.5 

.0457 

EJB 

+5 

7.5 

82.5 

.0303 

Similarly,  F was  calculated  for  p = 0.4,  0.2  and  0. 1 


S 

F 

.6 

0.0051 

.4 

0.0055 

2 

0.0130 

.1 

0.0238 

F = 0.0007  + .002326  p 


Steps  for  implementing  a wind  velocity  change  into  the  prediction  equation: 
1 . Calculate  F for  a given  p 
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1 


2.  Calculate  V£Q  = |(F)  (Vw)  + i .0]  [V  - Vw] 

where  = wind  component 

tail  wind  = -ve 
head  wind  = + ve 
V = baseline  V 

3-  V£q  ~y  ^or  Ui>lllS  the  prediction  equation 


105 


APPENDIX  B 

CALCULATION  OF  WEIGHT  CORRECTION  FACTOR 


Example:  (Variation  of  weight  alone): 


B-52,  p = 0.6,  Baseline  V = 90  Knots,  S = 1779\W  = 290,000  Lbs. 
Weight  only  variation  data  (simulator) 


\ 

W (lbs) 

AW  = 

W-W  (baseline) 

S (ft) 

(knots) 

(AW/AV) 

= (Aw/AVEQ  1 

Avg. 

(AW/AV) 

450,000 

160,000 

2822 

30.5 

5246 

370,000 

80,000 

2278 

15.2 

5263 

290,000 

0 

1779 

- 

- 

5150 

245,000 

- 45,000 

1517 

-8.7 

5172 

200,000 

- 90,000 

1246 

-18.3 

4918 

* The  AVgq  is  the  equivalent  velocity  change  from  baseline  (90  Knots)  needed  to  obtain 
the  distance,  S,  shown  in  the  previous  column,  and  is  obtained  from  the  velocity  only 
var.ation  data  shown  (table)  below. 


■6p 

Velocity  only  variation  data  (simulator) 

W 

V 

s 

Aveq. 

290,000 

90Kn 

1779 

0 

95 

1916 

5 

100 

2116 

10 

110 

2418 

20 

120 

2755 

30 

85 

1633 

-5 

80 

1457 

-10 

Similarly,  (AW/AV)  ( 4^}=  5640;  (AW/AV)  { 2(1)  = 5970: 
(AW/AV)  (!#l)  = 6088. 


X = #i 

Y = AW/AV 

.6 

5150 

.4 

5640 

.2 

5970 

.1 

6088 

Y = 6322-  1877  (X) 
or 


AV  = 


A0 

6322  - 1877  <p) 
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Steps  for  implementing  a weight  only  change  into  an  equivalent  velocity  only  change  for  use 
in  prediction  equation: 


Steps:  1 ) Calculate  (Aw/A v)  for  given  n. 

2)  Divide  total  Aw  by  (Aw/Av)  to  obtain  Av 

3)  Algebraically  add  Av  of  Step  2 to  baseline  v. 

\ 

I v = v + Av  for  higher  weights  and 

v = v - Av  for  lower  weights, 

4)  Then  use  (v  = v^asejjne  ± Av)  in  the  prediction  equation: 

1 : 

1 


APPENDIX  C 


CALCULATION  OF  PRE-BAKING  ROLL  DISTANCE  FOR  B-52  AIRPLANE 


A digital  computer  program  was  used  to  compute  the  distance  ‘S’  expressed  by  the  following 
integral 


THE  GROUND  RUN  OF  AN  AIRCRAFT  IS  GIVEN  BY: 
S = 2.852*  - (Vj  (V-Vw)dV 


j T Jv~  v^>dv 


S = distance  relative  to  the  ground  (ft.) 

W = weight  (lb) 

x = runway  gradient  (positive  upwards)  (radians) 

V = true  airspeed  (knots) 

Vy  = touch-down  speed  (knots) 

Vg  = brake  application  speed  (knots) 

D = total  drag  (lb)  = (aero  + drag  chute) 

G = (rolling)  friction  force  between  the  wheels  and  the  runway 
g ' = acceleration  due  to  gravity,  (32.2  ft/sec2) 

F = thrust  (lb) 

Vw  = Head  wind  component  (knots) 

*Tlus  is  the  conversion  factor  from  (knots)2  to  (ft/sec)2 
The  program  internally  computes  parameters  D,  G & F.; 


Example-  W = 450,000  lbs  CL=.305  Vw  = 0 

^rolling  = 02  x = 0 Vy  = 1 00  knots 

Cq  = .152  (Touchdown  to  Drag  Chute  deployment  speed  (135  knots) 
Cq  = .32 1 (135  knots  to  Vg  = 90  knots) 

160 

Then  IS]  = 3000  feet 

and  [ S)  ^ = 3375  feet 
160 

Total  |S)90  =6375  feet 

This  is  one  of  the  data  points  shown  in  figure  23, 


vw  = o 

Vy  = 100  knots 
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APPENDIX  D 

FORMULATION  OF  PREDICTION  MODEL 


I , DETERMINATION  OF  COMPONENT  EQUATIONS 

When  the  experimental  data  had  been  arranged  as  described  in  reference  (2).  Appendix  C, 
relationships  between  7r  j . the  term  containing  the  dependent  variable,  and  n 7.  ^3  and  ^4  in 
turn,  the  terms  with  independent  variables,  were  obtained  using  statistical  curve  fitting 
programs.  The  relationships  between  7r  j and  other  individual  7r  terms  are  called  component 
equations. 

Plots  were  prepared  of  7r  j versus  m,  7r  j versus  ^3,  and  7r  j versus  it 4 for  ah  atr;  lanes  using 
data  from  ASD-TR-77-6  Volume  if.  Section  XII.  An  example  of  these  plots  is  shown  (for 
F-i  1 1 data)  in  Figure  D-l,  This  helped  determine  the  general  form  ot  a relationship  that 
could  e'ost  between  rr  j and  m.  7T|  and  ^3,  and  so  on.  Figure  D-2  is  a flow  ciiait  depicting  the 
formulation  of  component  equations.  The  explanation  of  terms  linear  regression,  polynomial 
regression.- etc.  can  be  found  in  Appendix  D,  reference  (1). 

2.  DETERMINATION  OF  FUNCTIONS 

As  explained  in  Appendix  E,  reference  (1 ),  the  component  equations  were  combined  by  the 
product  method.  The  conditions  for  the  function  to  be  a product  and  the  equivalence  of 
calculated  values  of  C were  also  developed  in  the  referenced  Appendix,  The  same  procedure 
was  repeated  in  the  present  study.  Tables  D-l  through  D-4  show  detailed  calculations  for 
tests  of  validity. 

Having  established  the  type  of  component  equations,  a digital  computer  program  was  put 
together  to  systematically  calculate  n terms  (numerics),  formulate  component  equations, 
calculate  constant  and  calculate  correlation  errors. 
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Figure  D-2.— Equation  Flow  Chart 


Table  D- 1 .—Calculation  of  Validity  for  the  Function  to  be  a Product 


- 2 
(F) 

Fi 

b.9628 

2.4760 

6.1652 

2.4930 

29.4632 

5.4270 

60.3107 

7.9260 

1.3340 

1.1440 

2.8629 

1.693 

13.6678 

3.740 

66.4714 

8.035 

1 .9321 

1.411 

4.8753 

2.194 

27.6676 

5.122 

45.1315 

6.926 

Error 

Percentage 


1.410 

2.192 

5.283 

6.701 


-t •••••  '<'■■■  '4j».  *«•*••  » 


Table  D-3— Calculation  of  Validity  for  the  Function  to  be  a Product-Wet  Runway 


Airplane 

TT 

2 

-2 

(F) 

F1 

F2 

F x F, 
1 2 

(F)2 

Error 

Percentage 

B-52 

.23 

22.7529 

4.795 

4.788 

0.9910 

-0.89 

.175 

33.8491 

5.815 

5.669 

1 .0268 

2.68 

.127 

47.0870 

6.866 

6.492 

1 .0563 

5.63 

KC-135 

.166 

20.9123 

4.564 

4.515 

1.0148 

1.48 

.154 

25.0000 

5.008 

4.911 

1 .0165 

1.65 

.137 

32.6269 

5.712 

5.629 

1.0147 

1.47 

F-lll 

.231 

19.3864 

4.447 

4.375 

0.9964 

-0.36 

.218 

21 .8930 

4.694 

4.663 

1 .0002 

0.02 

.193 

28.2705 

5.352 

5.358 

0.9859 

-1.41 

I/C 

= (F) 

2 - [f  (S 

2-  v V] 

2 

TEST  OF  VALIDITY: 

(F-i)  - F 1T^»  ^4) 

j . (Fl) 

(F2)  = F (tt2,  Vy  tt4)  - 2 
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APPENDIX  E 


CALCULATION  OF  DENSITY  CORRECTION  FACTOR: 

Definitions. 

Standard  day  = 59°F,  sea  level,  p = .00238  lb-scc2/ft4 
Non-st'andard  day  = Hot  day  or  cold  day  as  follows: 

Hot  day  = 83°F/28°C,  p = .00189  ib-sec2/ft4 

Cold  day  = -60°F/-5 1°C,  p = .00309  lb-sec2/ft4 

Example:  KC-135;0.!p 

Component  Equation  25  (unmodified)  had  the  form:; 

(jt,  ) = 82.268669  (ir4)  --200863 

When  this  equation  was  used  to  predict  stopping  distances  for  hot  and  cold  day  cases  the 
correlation  error  exceeded  the  desired  ± 5%  (see  the  following  tabular  data). 


Actual 

Predicted 

% 

Condition 

fir4) 

(Wj) 

fir,) 

Error 

Hot  day 

mm 

9.131 

8.453 

7.420 

Cold  day 

HI 

6.976 

7.658 

-9.776 

A modified  exponent  was  calculated  for  the  term  (tt4)  such  that  the  correlation  error  would 
be  minimal. 

Let  modified  exponent  = B 

/ Actual  (ffj  X/ 

Then  for  no  error ; B = log  ^ 7’“  hog  (ir4) 

i.e,  for  a hot  day; 

bHOT  = 108  (9- 1 3 1 /8  2.2687)/log  (855 1 7)  = -.  1 93574 
A exponent  = .200863  -.193574=  .007289 
Similarly  = -.208263 

and  AeXp0  = .007400 
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This  Aexpo  cannot  be  directly  added  to  the  7r4  exponent  since  the  ff4  term  is  also  used  to 
predict  affects  of  varying  V and  Fe.  i.e.  an  explicit  form  of  Aexp0  for  density  change 
A p should  be  calculated.  This  is  illustrated  below: 


Condition 

P 

(density) 

hp 

Exponent 

^x/Ap 

(avg.) 

Std  day  ' 

.00238 

- 

- 

- 

- 

Hot  day 

.00189 

.00049 

.007289 

14.875 

12.65 

Cold  day 

.00309 

.00071 

.007400 

10.423 

The  modified  component  equation  is:. 


j : 


I2.65A p -.200863 

<7T|)  = 82.268669  (jt4)l 


whereAp  = pstd-pnon-Std 


(25) 


I 

i 

i 


4 
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